NASA TECHNICAL NOTE 


NASA TN D-8430 



CO 

CO 


I 





i N 4S k N -L*- b 4 3 0 ) hEVISiiD FQBTaAN FBCGbAfl 

ZQ& CALCULATING VbLOCxlx.cS AND STB cAMLlNbS 
JN T Hb hUb-SHBOUL MIDCtiA NNbl STcLaM SUBFACb 

CF AN AXIAL-/ BAL1AL-, 0* fllXBL-FLGW 

1 UBbOdAC alisb CE ANNULAE DUCT. 1* UbFB 5 ril/01 


N77- 18992 

Ht AOU 
MF AO/ 

Uaclas 

21299 



REVISED FORTRAN PROGRAM FOR 
CALCULATING VELOCITIES AND STREAMLINES 
ON THE HUB-SHROUD MIDCHANNEL STREAM 
SURFACE OF AN AXIAL-, RADIAL-, OR 
MIXED-FLOW TURBOMACHINE OR ANNULAR DUCT 


I - User’s Manual 


Theodore Katsanis and William D. McNally 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS ANO SPACE ADMINISTRATION • WASHINGTON, 0. C. - MARCH 1977 




1. Report No. 2. Government Accession No 

NASA TN D-8430 

3. Recipient's Catalog No. 

4. Title and Subtitle REVISED FORTRAN PROGRAM FOR CALCULATING 

VELOCITIES AND STREAMLINES ON THE HUB-SHROUD MIDCHANNEL 
STREAM SURFACE OF AN AXIAL-, RADIAL-, OR MIXED- FLOW 
TURBOMACHINE OR ANNULAR DUCT 
I - USER’S MANUAL 

5. Report Date 

March 1977 

6. Performing Organization Code 

7. Author(s) 

Theodore Katsanis and William D. McNally 

8. Performing Organization Report No, 

E -8734 

10. Work Unit No. 

505-04 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Supersedes NASA TN D-7343 (FORTRAN Program for Calculating Velocities and Streamlines 
on the Hub-Shroud Mid-Channel Flow Surface of an Axial- or Mixed-Flow Turbomachine, 1973) 


16. Abstract 

A FORTRAN -IV computer program has been developed that obtains a detailed subsonic or shock - 
free transonic flow solution on the hub -shroud midchannel stream surface of a turbomachine. 

The blade row may be fixed or rotating, and the blades may be twisted and leaned. Flow may be 
axial, mixed, or radial. This program is a revision of a previous program and this report 
supersedes NASA TN D-7343. Upstream and downstream flow variables may vary from hub to 
shroud, and provision is made to correct for loss of stagnation pressure. The results include 
velocities, streamlines, and flow angles on the stream surface as well as approximate blade 
surface velocities. Subsonic solutions are obtained by a finite -difference, stream -function so- 
lution. Transonic solutions are obtained by a velocity -gradient method that uses information 
from a finite -difference, stream -function solution at a reduced mass flow. 


17. Key Words {Suggested by Author(sl) 

Meridional plane; Turbomachine flow; Midchannel 
stream surface; Axial -flow turbomachine; Mixed - 
flow turbomachine; Transonic flow; Radial -inflow 
turbine; CentrifugaL compressor 


18. Distribution Statement 

Unclassified - unlimited 
STAR Category 01 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif, {of this page) 

Unclassified 


21. 


No. of Pages 

105 


22, Price* 

A06 


* For sale by the National Technical Information Service. Springfield. Virginia 22161 













CONTENTS 


Page 

SUMMARY 1 


INTRODUCTION 2 


METHOD OF ANALYSIS 5 

BASIC ASSUMPTIONS 5 

SOLUTION BY COMBINATION OF METHODS 6 

SUBSONIC STREAM- FUNCTION SOLUTION 6 

TRANSONIC VELOCITY-GRADIENT APPROXIMATE SOLUTION 8 

BLADE SURFACE VELOCITIES 9 

APPLICATION OF PROGRAM 9 


DESCRIPTION OF INPUT AND OUTPUT 10 

INPUT 11 

Input Dictionary 17 

Special Instructions for Preparing Input 25 

(a) Units of measurement 26 

(b) Damping factors FNEW and DNEW 26 

(c) Hub and shroud flow channel geometry 26 

(d) Orthogonal mesh 27 

(e) Upstream and downstream flow conditions 27 

(f) Mean blade surface and thickness coordinates 28 

(g) How to specify points for spline curves 29 

(h) Requests for output data 29 

(i) Incompressible flow 31 

(j) Straight infinite cascade 31 

(k) Choosing a value for REDFAC 31 

( l ) Choosing a value for ANGROT 31 

OUTPUT 31 

Printed Output 32 

Plotted Output 52 

Error Messages 6€ 

NUMERICAL EXAMPLE 70 


iii 


PRECEDING PAGE BLANK NOT FILMED 



APPENDIXES 

A - GOVERNING EQUATIONS 73 

B - DERIVATION OF STREAM- FUNCTION EQUATION 77 

C - DERIVATION OF VELOCITY- GRADIENT EQUATION 85 

D - LOSS CORRECTIONS 90 

E - DEFINING REDUCED- MASS- FLOW PROBLEM 92 

F - INCIDENCE AND DEVIATION CORRECTIONS 93 

G - BLADE SURFACE VELOCITIES AND BLADE- TO-BLADE 

AVERAGE DENSITIES 95 

H - SYMBOLS 98 

REFERENCES 101 



REVISED FORTRAN PROGRAM FOR CALCULATING VELOCITIES AND 
STREAMLINES ON THE HUB-SHROUD MIDCHANNEL STREAM 
SURFACE OF AN AXIAL-, RADIAL-, OR MIXED-FLOW 
TURBOMACHINE OR ANNULAR DUCT 
I - USER'S MANUAL* 

by Theodore Katsanis and William D. McNally 
Lewis Research Center 

SUMMARY 

A FORTRAN-IV computer program has been developed that obtains a detailed sub- 
sonic or transonic flow solution on the hub-shroud midchannel stream surface of a 
single blade row of a turbomachine. A solution can also be obtained for an annular 
duct without blades. The flow must be essentially subsonic, but there may be locally 
supersonic flow. The solution is for two-dimensional, adiabatic shock-free flow. The 
blade row may be fixed or rotating, and the blades may be twisted and leaned. The 
flow may be axial, mixed, or radial. Upstream and downstream flow conditions can 
vary from hub to shroud, and provision is made for an approximate correction for loss 
of stagnation pressure. Viscous forces are neglected along solution mesh lines running 
from hub to tip. 

The present program is a revision of a previous program and this report super- 
sedes NASA TN D-7343. The primary revisions are to extend the program to handle 
nonaxial flows without restriction, to handle annular ducts without blades, to allow for 
any specified streamwise loss distribution, and to make numerous detailed improve- 
ments for more accurate and efficient calculations. 

The basic analysis is based on the stream function and consists of the solution of 
the simultaneous, nonlinear, finite- difference equations of the stream function. This 
basic solution, however, is limited to strictly subsonic flow. When there is locally 
supersonic flow, a transonic solution must be obtained. The transonic solution is ob- 
tained by a combination of a finite- difference, stream-function solution and a velocity- 
gradient solution. The finite- difference solution at a reduced mass flow provides infor- 
mation that is used to obtain a velocity- gradient solution at the full mass flow. 

The program is reported in two volumes, with part I as the user's manual and 
part II as the programmer's manual. This report, part I, contains all the information 
necessary to use the program as it is. It explains the equations involved and the method 
of solution and gives a numerical example to illustrate the use of the program. Part n 
includes the complete program listing and a detailed program procedure. 

*Supersedes NASA TN D-7343. 



INTRODUCTION 


The design of blades for compressors and turbines ideally requires methods for 
analyzing unsteady, three-dimensional, turbulent viscous flow through a turbomachine. 
Clearly, such solutions are impossible at the present time, even on the largest and fast- 
est computers. The usual approach at present is to analyze only steady flows and to 
separate inviscid solutions from viscous solutions. Three-dimensional inviscid solu- 
tions are just beginning to be obtained with the present generation of computers, but they 
require excessive computer time. So, at present, inviscid analyses usually involve a 
combination of several two-dimensional solutions on intersecting families of stream 
surfaces to obtain what is called a quasi- three-dimensional solution. 

Since there are several choices of two-dimensional surfaces to analyze and many 
ways of combining them, there are many approaches to obtaining a quasi-three- 
dimensional solution. Most two-dimensional solutions are either on a blade-to-blade 
surface of revolution (Wu's S^ surface, ref. 1) or on the meridional or midchannel 
stream surface between two blades (Wu's Sg surface). However, when three- 
dimensional effects are most important, significant information can often be obtained 
from a solution on a passage cross-sectional surface (normal to the flow). This is 
called a channel solution (fig. 1) . 

In this report a solution to the equations of flow on the meridional surface is 
carried out. This solution surface is chosen when the turbomachine under consideration 
has significant variation in flow properties in the hub-shroud direction, especially when 
input is needed for use in blade-to-blade calculations. The solution can be obtained 
either by the qua si- orthogonal method, which solves the velocity- gradient equation from 
hub to shroud on the meridional stream surface (ref. 2) , or by a finite-difference me- 
thod, which solves a finite- difference equation for stream function on the same stream 
surface. The quasi-orthogonal method is efficient in many cases and can obtain solu- 
tions into the transonic regime. However, there is difficulty in obtaining a solution 
when blade aspect ratios are above 1. Difficulties are also encountered with curved 
passages and low-hub- tip- ratio blades. For such cases, the most promising method is 
the finite-difference solution, but this solution is limited to completely subsonic flows. 

Finite-difference programs for flow on the midchannel surface of a turbomachine 
have been reported in the literature. However, many of these programs are proprie- 
tary or are of limited generality. The program reported herein is very general and has 
been thoroughly tested and refined as the result of extensive usage at the NASA Lewis 
Research Center. 

The program described in this report uses both the finite- difference and the quasi- 
orthogonal (velocity gradient) methods, combined in a way that takes maximum advan- 
tage of both. The finite- difference method is used to obtain a subsonic-flow solution. 
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Figure 1. -Two-dimensional analysis surfaces in a turbomachine. 


The velocity-gradient method is then used if necessary to extend the range of solutions 
into the transonic regime. 

A computer program called MERIDL has been written to perform these calcula- 
tions. This program was written for an axial-, mixed-, or radial-flow turbomachine 
blade row, either a compressor or turbine, or for an annular duct. Upstream and down- 
stream flow conditions can vary from hub to shroud. The solution is for compressible, 
shock-free flow or for incompressible flow. An approximate correction for loss of 
stagnation pressure through the blade row is provided. The blade row may be either 
fixed or rotating. The blades may be twisted and leaned and can have high aspect ratio 
and arbitrary thickness distribution. 

The solution obtained by this program also provides the information necessary for 
a more detailed blade-shape analysis on blade- to-blade surfaces (fig. 1) . A useful pro- 
gram for this purpose is TSONIC (ref. 3) . Information needed to prepare all the input 
for TSONIC is calculated and printed by MERIDL. 

The MERIDL program reported herein is a revision of the program described in 
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references 4 and 5. Two types of changes were made: first, extension of the program's 
capability to handle cases beyond those originally offered; and second, revisions to im- 
prove the accuracy and efficiency of the program. Although the input form has been ex- 
tended to handle additional input where required, any input that was satisfactory for the 
original MEREDL program is still satisfactory for the revised MERIDL program. The 
following major extensions and revisions were made to the program (additional internal 
changes are documented in part II): 

(1) The program has been extended to handle nonaxial flow configurations, including 
both radial turbines and centrifugal compressors. 

(2) The program has been extended to handle an annular duct without blades. 

(3) The program has been extended to permit the user to specify an arbitrary 
streamwise distribution of loss within the blades. This is in addition to the original 
provision for hub-shroud loss distribution. 

(4) The program has been modified so that the blade thickness can be specified pre- 
cisely by a set of tangential (circumferential) thickness coordinates. The original pro- 
gram required specification of thickness normal to the mean camber line on an input 
blade section. This normal thickness was influenced by blade lean, camber, and non- 
parallel blade surfaces and was difficult to specify accurately for some blade shapes. 

(5) If desired the leading- and trailing-edge, mean-camber- line tangent angles can 
be specified as input. This simplifies the specification of some blade shapes. 

(6) Output quantities have been added to station- line output to give absolute velocity 
components and to give static as well as absolute and relative total temperature, den- 
sity, and pressure . 

(7) Several informational messages have been added to the output. 

(8) Additional error messages have been provided. 

(9) Upstream and downstream boundary conditions have been changed to improve 
both the convergence and the quality of the solution near these boundaries. 

(10) Interpolation and calculation procedures near the leading and trailing edges have 
been improved to give better convergence and smoother solutions in these regions. 

(11) Numerous small changes have been made to improve the accuracy and efficiency 
of the program. 

The MERIDL program has been implemented on the NASA Lewis time-sharing 
IBM-TSS/360-67 computer. For the numerical example of this report, the storage of 
variables required 60 000 words for a 21 x 41 grid of 861 points. The amount of varia- 
ble storage could be easily reduced by equivalencing of variables or by using a coarser 
mesh. Storage requirements for the program code depend on the computer system and 
compiler being used. Run times for the program range from 3 to 15 minutes on 
EBM-TSS/360-67 equipment. 

The MERIDL program is reported in two volumes, with part I as the user's manual 
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and part II (ref. 6) as the programmer's manual. This report, parti, contains all the 
information necessary to use the program. It explains the method of solution, de- 
scribes the input and output, gives a numerical example to illustrate the use of the pro- 
gram, and derives the equations used (in appendixes A to G.) Symbols are defined in 
appendix H. Part II includes a complete program listing, detailed program procedure, 
and appendixes that derive the special numerical techniques used. 


METHOD OF ANALYSIS 
BASIC ASSUMPTIONS 

It is desired to determine the flow distribution through a stationary or rotating 
cascade of blades on a midchannel hub-shroud stream surface. The following simpli- 
fying assumptions are used in deriving the equations and in obtaining a solution: 

(1) The flow relative to the blade is steady. 

(2) The fluid is a perfect gas with constant specific heat C . 

(3) The only forces along a hub- shroud orthogonal mesh line are those due to 
momentum and pressure gradient. 

(4) There is no heat transfer. 

(5) The midchannel surface is a stream surface that has the same shape as the 
blade mean-camber surface, except near the leading and trailing edges, where an arbi- 
trary correction is made to match the free-stream flow. 

(6) The velocity varies linearly between blade surfaces. 

(7) The relative stagnation pressure loss is known through the blade row. 

(8) The upstream and downstream boundaries of the solution region are orthogonal 
to the streamlines. 

The flow may be axial, mixed, or radial. Whirl, stagnation pressure, and stagna- 
tion temperature may vary from hub to shroud, both upstream and downstream of the 
blade row. The blade row may be either fixed or rotating, with leaned and twisted 
blades; or there may be no blades at all. Within the given assumptions, no terms are 
omitted from the equations . 

In connection with assumption 3, the viscous forces along a hub-shroud orthogonal 
mesh line are neglected, since these forces are usually very small. The viscous for- 
ces in the streamwise direction are much larger. The effects of these streamwise 
viscous forces are considered indirectly by specifying a streamwise total-pressure- 
loss distribution. 
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SOLUTION BY COMBINATION OF METHODS 


A flow analysis on the meridional stream surface can be obtained either by the 
velocity-gradient method or by the finite-difference method. The finite- difference 
method is limited to subsonic flow; the velocity- gradient method is limited to relatively 
low-aspect- ratio blades. The most accurate solution is obtained by the finite- 
difference technique, so that this method is used where possible (i. e. , for subsonic 
flow). With locally supersonic flow, the finite- difference solution is first obtained at 
a reduced mass flow (appendix E) for which the flow field is completely subsonic. The 
streamline curvatures and flow angles throughout the passage that are obtained from 
this solution provide the information necessary to obtain an approximate velocity- 
gradient solution at full mass flow, regardless of aspect ratio. 


SUBSONIC STREAM- FUNCTION SOLUTION 

The stream- function equation (eq. (Al), appendix A) is a partial differential equa- 
tion on a midchannel hub-shroud stream surface (assumption 5, p. 5). It is in one un- 
known (the stream function) as a function of two variables and is derived from Wu's 
momentum equation (eq. (96), ref. 1) on what he calls an S 2 surface. Equation (A 1) is 



Section A-A 
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nonlinear but can be solved iteratively by the finite- difference method when the flow is 
completely subsonic. 

A finite region (as indicated in fig. 2) is considered for the solution of equa- 
tion (Al) . It is assumed that the upstream and downstream boundaries are sufficiently 
far from the blade so as to have a negligible effect on the solution. Equation (Al) is 
elliptic for subsonic flow. Therefore, when the flow is entirely subsonic, equation (Al) 
can be solved when proper boundary conditions are specified on the entire boundary of 
the region. These conditions are the values of the stream function or its normal de- 
rivative on all four boundaries. The stream function is zero on the hub and 1 at the 
shroud. On the upstream and downstream boundaries, the derivative of the stream 
function normal to the boundary is assumed to be zero. This is equivalent to assuming 
that the upstream and downstream boundaries are orthogonal to the streamlines (as- 
sumption 8, p. 5). 

The numerical solution of equation (Al) is obtained by the finite— difference method. 
A grid must be used for the finite- difference equations. The type of grid that was 
chosen is an orthogonal mesh that is generated by the program. The method used to 
generate the grid follows that reported in reference 7. The space between the hub and 
the shroud is divided into equal increments along several hub- shroud lines. Spline 
curves are then fit through the resulting points to obtain the streamwise orthogonals 
(fig. 3). The normal orthogonals are obtained by a predictor- corrector technique. 

This technique is analogous to the second-order Runge-Kutta method for solving ordi- 
nary differential equations, which is known also as the improved Euler method or 
Heun's method (ref. 8). The distances along orthogonal mesh lines are s and t. The 
s— distance is in the streamwise direction and the t— distance is normal to this, as indi- 
cated in figure 3. With the mesh determined, the finite-difference equations on the 
orthogonal mesh are given in part II of this report (ref. 6). Note that the orthogonal 
mesh is not used as a coordinate system. 



Figure 3. - Orthogonal finite-difference mesh on solution region. 
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The finite- difference equations are nonlinear since the original equation (Al) is non- 
linear. These equations can be solved iteratively. On the first iteration an initial den- 
sity is assumed; this linearizes some of the terms. The remaining nonlinear terms are 
omitted for the first iteration so that the finite- difference equations are entirely linear- 
ized. These linearized equations are then solved to obtain the first approximate solu- 
tion for the stream function. This solution provides information that is used to obtain a 
better estimate of the density and an estimate of the other nonlinear terms. The equa- 
tions are then solved again to obtain an improved solution. This process is repeated, 
and by iteration a final converged solution can be obtained if the flow is subsonic. 

For each step of this iteration the linearized finite-difference equations must be 
solved. The method used to solve the equations is successive overrelaxation (ref. 9) 
with an optimum overrelaxation factor. Since this is also an iterative method, we have 
two levels of iteration. The overrelaxation is performed in the "inner iteration," and 
the corrections to the nonlinear terms are made in the "outer iteration." The inner 
iteration is internal to the program, so the only iteration apparent to the user is the 
outer iteration. When an iteration is mentioned in this report, it is always the outer 
iteration. 

After the stream function is obtained, the velocity distribution is found by numeri- 
cal partial differentiation of the stream function and by using equations (A5) and (A6). 

The details of the numerical procedure and programming technique are described in 
part II (ref. 6) . 


TRANSONIC VELOCITY- GRADIENT APPROXIMATE SOLUTION 


For the case where there is locally supersonic flow, equation (Al) is no longer 
elliptic in the entire region but is hyperbolic in the region of supersonic flow (ref. 10). 
This changes the boundary conditions and means that there will probably be shock 
losses in going from supersonic to subsonic flow. The finite-difference method cannot 
be used with locally supersonic flow. However, an approximate solution can be ob- 
tained by getting a reduced-flow solution with the finite-difference method and extending 
this to the full flow by using the velocity- gradient method. This technique is described 
in reference 3. 

The velocity- gradient equations are equations (A 7) to (All). Equation (A7) is 
solved as an initial-value problem, where the velocity W is specified at the hub for 
any given vertical mesh line running from hub to tip. By finding several solutions for 
varying values of W at the hub, a solution satisfying the specified mass flow 
(eq. (A 12)) will be found. When equation (A7) has been solved, subject to giving the 
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correct mass flow, for every hub-shroud mesh line in the region, the entire velocity 
distribution at full mass flow has been obtained. 


BLADE SURFACE VELOCITIES 

The solution that is obtained by either the finite- difference or velocity- gradient 
method is for the midchannel surface between the blades. When there are blades, the 
blade surface velocities are of greater interest. These can be estimated since the 
blade loading is dependent on the rate of change of whirl. By assuming a linear varia- 
tion of velocity between blade surfaces, equation (A13) can be derived for calculating 
the blade surface velocities. 


APPLICATION OF PROGRAM 


The program can be used both for analysis and as a design tool. When it is used 
for design, other programs should be used with this program. For axial compressors, 
reference 11 describes a program that will give blade mean- camber-line coordinates 
and thicknesses for an axial compressor blade. This blade design can be checked by 
using the MERIDL program to analyze the flow distribution in detail. Usually, changes 
must be made to the blade design to achieve a desirable flow distribution. These 
changes may involve more than just the blade shape; for example, hub and shroud pro- 
file, inlet and outlet whirl distribution, and loss distribution may have to be changed. 

Of course, the accuracy of the MERIDL solution depends on the accuracy of the bound- 
ary conditions used. 

When a reasonable flow pattern is achieved by the MERIDL flow analysis on the 
midchannel stream surface, more-detailed blade surface velocities can be obtained by 
flow analyses on various blade-to-blade stream surfaces. A useful program for this 
purpose is TSONIC (ref. 3). Most of the information required to compute input for 
TSONIC is calculated and printed directly by MERIDL. Item 12 of the section Printed 
Output (p. 50) explains how to obtain and use this output from MERIDL. Further 
changes in blade shape or whirl distribution may be considered at this time. Refer- 
ence 12 (ch. VTI) gives information on incidence and deviation for good designs and for 
off-design conditions. 

For cases where the flow is well guided in the channel but has large variations, 
both from blade to blade and from hub to shroud, the CHANEL program (ref. 13) is 
useful. The CHANEL program obtains a solution on a channel cross-section surface. 
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The CHANEL program is particularly useful for calculating choking mass flow through 
a blade row. 


DESCRIPTION OF INPUT AND OUTPUT 


The principal block of input required for the program is a geometric description of 
the blade row to be analyzed. The geometry is given in cylindrical coordinates, as 
shown in figure 4. The blade shape is given as a series of blade sections from hub to 
shroud. Each blade section is described by a set of z and r mean- camber-line 
coordinates and two sets of blade- shape coordinates on any general, smooth surface of 



Figure 4. - Cylindrical coordinate system and velocity components. 
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revolution. Other inputs include upstream and downstream flow conditions, a geo- 
metric description of hub and shroud, appropriate gas constants, operating conditions 
such as mass flow and rotational speed, and a description of the finite- difference solu- 
tion mesh. There are also input variables to indicate the desired blocks of output data. 

Ouptut is given at any or all of three principal geometric locations: (1) on all mesh 
points of the orthogonal solution mesh, (2) along user- designated streamlines through 
the blade row, or (3) along user- designated lines from hub to shroud (station lines). 
Output consists mainly of z- and r- coordinates and velocity components and angles. 

At station lines more output is given, including density, temperature, and pressure. 


INPUT 

Figure 5 shows the input variables as they are punched on the data cards. The first 
input data card is for a title, which serves for problem identification. Any information 
may be put in the 80 columns of this card. 

All the numbers on the three input cards beginning with MBI, LSFR, and IMESH are 
integers (no decimal point) in a five- column field (fig. 5). These must be all right ad- 
justed. The input variables on all other data cards are real numbers (punch decimal 
point) in 10- column fields. 

Figure 5 indicates that several options exist for specifying upstream and down- 
stream flow conditions. First, the user can specify these conditions either as a func- 
tion of the stream function (SFIN and SFOUT) or as a function of radius (RADIN and 
RADOUT). Next, the user can specify either whirl (LAMIN and LA.MOUT) or absolute 
tangential velocity (VTHIN and VTHOUT) . At the downstream input station either abso- 
lute total pressure (PROP) or fractional loss of total pressure (LOSOUT) may be given. 
Finally, a loss distribution (PERLOS) within the blade may be specified. If the loss 
distribution through the blade is not given, the loss is assumed to vary linearly between 
the leading and trailing edges . 

Input variables are both geometric and nongeometric. The geometric input varia- 
bles are shown in figures 6 to 10. Further information concerning the input variables 
is given in the section Special Instructions for Preparing Input (pp. 25 to 31) . 
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Figure 5. - Input form. (Cards denoted by asterisks are optional. See Input Dictionary. ) 















Input Dictionary 


The input variables are described in terms of a consistent set of SI units: newtons, 
kilograms, meters, joules, kelvins, and seconds. The program, however, will run 
with input in any consistent set of units. 

The input variables, in the order that they appear in figure 5, are the following: 


GAM 

AR 

MSFL 

OMEGA 

REDFAC 


VELTOL 


FNEW 


DNEW 


MBI 


Specific-heat ratio, y 

Gas constant, J/ (kg) (K) 

Total mass flow through entire circumferential annulus of blade 
row, kg/ sec 

Rotational speed, u>, rad/sec. Note that cu is negative if rotation 
is in opposite direction of that shown in fig. 4. 

Factor by which mass flow (MSFL) must be reduced in order to 
assure subsonic flow throughout flow passage. REDFAC may be 
left blank, in which case a value of 1.0 will be used. See sec- 
tion (k), p. 31. 

Convergence tolerance on maximum velocity change in each iter- 
ation, over all mesh points, for reduced- mass- flow solution. 
VELTOL may be left blank, and a value of 0.01 will be used by 
the program. Whatever value is given is multiplied by the mini- 
mum of FNEW and DNEW before it is printed and subsequently 
used by the program. A value of 0.001 for VELTOL is a tight 
tolerance, 0.01 is a medium tolerance, and 0.1 is a loose toler- 
ance. 

Damping factor on calculation of F^. from iteration to iteration. 

A value of 0. 5 is suggested for FNEW. FNEW may be left blank, 
in which case the program will use a value of 0. 5. See sec- 
tion (b), p. 26. 

Damping factor on calculation of O(rV 0 )/3t within blade row from 
iteration to iteration. DNEW may be left blank, and the program 
will use a value of 0. 5. DNEW does not have to be equal to 
FNEW. See section (b), p. 26. 

Number of vertical mesh lines from left boundary of orthogonal 
mesh (ZOMIN) to point of first mesh-size change (ZOMBI). Use 
MBI = 0 if there are no blades. See fig. 8 and section (d), p. 27. 
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MBO 

MM 

MHT 
NBL 
Nil UB 
NTIP 
NIN 

NOUT 

NBLPL 

NPPP 

NOS TAT 


NSL 


Total number of vertical mesh lines from left boundary of orthogo- 
nal mesh (ZOMIN) to point of second mesh-size change (ZOMBO). 
Use MBO = 0 if there are no blades. See fig. 8 and section (d), 
p. 27. 

Total number of vertical mesh lines from left to right boundaries 
of orthogonal mesh (ZOMIN to ZOMOUT), maximum of 100. See 
fig. 8 and section (d), p. 27. 

Total number of horizontal mesh spaces from hub to shroud of or- 
thogonal mesh, maximum of 100. See fig. 8 and section (d), p. 27. 

Number of blades in total circumference of blade row. Use 
NBL = 1 if there are no blades. 

Number of spline points given in ZHUB and RHUB arrays, maxi- 
mum of 50. See fig. 6 and section (c), p. 26. 

Number of spline points given in ZTIP and RTIP arrays, maximum 
of 50. See fig. 6 and section (c), p. 26. 

Number of data points given in upstream arrays of flow properties 
(SFIN, RADIN, TIP, PRIP, LA MIN, VTHIN), maximum of 50. 

See fig. 7 and section (e), p. 27. 

Number of data points given in downstream arrays of flow proper- 
ties (SFOUT, RADOUT, PROP, LOSOUT, LAMOUT, VTHOUT), 
maximum of 50. See fig. 7 and section (e), p. 27. 

Number of blade planes or blade sections on which data (ZBL, 

RBL, etc.) are given to describe mean blade shape and blade 
thickness, maximum of 50. See fig. 6 and section (f) , p. 28. 

Number of data points per blade section or blade plane in ZBL, 
RBL, etc. , arrays, maximum of 50. See fig. 10 and sec- 
tion (f), p. 28. 

Number of stations from hub to shroud (located by coordinates in 
ZHST and ZTST) at which output is desired, maximum of 50. 

See fig. 9 and section (h), p. 29. NOSTAT may be left blank, in 
which case no input cards should be included for ZHST, ZTST, 
RHST, and RTST arrays. 

Number of streamlines from hub to shroud (designated by values 
in FLFR) at which output is desired, maximum of 50. See fig. 9. 
NSL may be left blank, in which case no cards should be included 
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NLOSS 


LSFR 

LTPL 

LAMVT 

LROT 

LB LAD 


LETEAN 

ANGROT 

ZOMIN 

ZOMBI 


for FLFR array. If NSL is left blank, theprogram will set it 
equal to 11 and print requested streamline output on 11 stream- 
lines that vary by 10 percent of total flow (i.e. , 0, 10, 20 

100 percent) . 

Number of points at which fractional loss (PERLOS) is specified, 
maximum of 50. If NLOSS is zero or left blank, the program will 
use a linear loss distribution in the streamwise direction. In this 
case, no input cards should be included for PERCRD or PERLOS 
arrays. 

Integer (0 or 1) indicating whether upstream and downstream flow con- 
ditions are given as a function of stream function (0) or radius (1) . 

Integer (0 or 1) indicating whether downstream total pressure (0) 
or fractional loss of stagnation pressure (1) is given as input. 

Integer (0 or 1) indicating whether upstream and downstream whirl 
(0) or tangential velocity (1) is given as input. 

Integer (0 or 1) indicating whether coordinate rotation (ANGROT) 
is necessary (1) for hub and shroud spline fit curves. See sec- 
tion ( l ), p. 31. 

Integer (0, 1, or 2) indicating which two blade- shape coordinates 
are given as input. If LB LAD = 0, mean blade shape 8 (THBL) 
and normal thickness t (TNBL) are given. If LB LAD = 1, 
mean blade shape 9 (THBL) and tangential 0-thickness tV r 
(TTBL) are given. If LBLAD = 2, upper blade surface 8 ^ 

(TH1BL) and lower blade surface 8 ^ (TH2BL) are given. See 
fig. 10 and section (f), p. 28. 

Integer (0 or 1) indicating whethe r leading- and trailing- edge mean 
blade- shape angles and /3 (BETALE and BETATE) are 
specified (1). See fig. 10 and section (f), p. 28. 

Rotation angle of axis used for all streamwise spline fit curves, 
deg. Omit this card if LROT is (0). See section (l) , p. 31. 

z- Coordinate of intersection of left boundary of orthogonal mesh 
with hub profile, m. See figs. 6 and 8 and section (d), p. 27. 

z- Coordinate of intersection of vertical mesh line with hub profile 
where first change in mesh spacing occurs (MBI), m. See figs. 

6 to 8 and section (d), p. 27. Set ZOMBI and ZOMBO to zero 
for no blades case. 
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ZOMBO 


ZOMOUT 

ROM IN, ROMBI, 
ROMBO, and 
ROMOUT 

ZIIUB 

RIIUB 

ZTIP 

RTIP 

ZHIN 

ZTIN 

RHIN. RTIN 
SPIN 
RADIN 
TIP 


z- Coordinate of intersection of vertical mesh line with hub pro- 
file where second change in mesh spacing occurs (MBO), m. 

See figs. 6 and 8 and section (d), p. 27. 

z- Coordinate of intersection of right boundary of orthogonal mesh 
(MM) with hub profile, m. See figs. 6 and 8 and section (d), p. 27. 

r-Coordinates corresponding to ZOMIN, ZOMBI, ZOMBO, and 
ZOMOUT, m. Leave these spaces blank if LROT = 0. 

Array of z- coordinates of input points defining hub or bottom 
boundary of flow channel, m. See fig. 6 and section (c) , p. 26. 

Array of r-coordinates of input points defining hub or bottom 
boundary of flow channel, m. See fig. 6 and section (c) , p. 26. 

Array of z- coordinates of input points defining shroud or top 
boundary of flow channel, m. See fig. 6 and section (c), p. 26. 

Array of r-coordinates of input points defining shroud or top 
boundary of flow channel, m. See fig. 6 and section (c), p. 26. 

z-Coordinate of intersection with hub profile of line on which up- 
stream flow conditions are given, m. Leave this entire card 
blank if LSFR = 0, LAMVT = 0, and MBI ^ 0. See fig. 7 and 
section (e), p. 27. 

z-Coordinate of intersection with shroud profile of line on which 
upstream flow conditions are given, m. See fig. 7 and sec- 
tion (e), p. 27. 

r-Coordinates corresponding to ZHIN and ZTIN, m. Leave these 
spaces blank if LROT = 0. If RADIN is given as input 
(LSFR = 1) , RHIN and RTIN cannot be equal to each other. 

Array of values of stream function for input points from hub to 
shroud along line on which upstream flow conditions are given. 
SFIN is given when LSFR = 0. See fig. 7 and section (e), p. 27. 

Array of r-coordinates of input points along line from hub to 
shroud on which upstream flow conditions are given, m. RADIN 
is given when LSFR = 1. See fig. 7 and section (e), p. 27. 

Array of absolute total temperatures T! at input points along line 
from hub to shroud on which upstream flow conditions are 
given, K. See fig. 7 and section (e), p. 27 
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PRIP 


LAMIN 


VTHIN 


ZHOUT 


ZTOUT 


Array of absolute total pressures p! at input points along line 
from hub to shroud on which upstream flow conditions are given, 
N/m 2 . See fig. 7 and section (e), p. 27. 

Array of values of absolute whirl ^rV^ at input points along line 

from hub to shroud on which upstream flow conditions are given, 
m 2 /sec. LAMIN is given when LAMVT= 0. See fig. 7 and sec- 
tion (e), p. 27. 

Array of values of absolute tangential velocity at input 

points along line from hub to shroud on which upstream flow con- 
ditions are given, m/ sec. VTHIN is given when LAMVT = 1. 

See fig. 7 and section (e), p. 27. 

z-Coordinate of intersection with hub profile of line on which 
downstream flow conditions are given, m. Leave this entire 
card blank if LSFR = 0, LAMVT = 0, and MBI / 0. See fig. 7 
and section (e), p. 27. 

z-Coordinate of intersection with shroud profile of line on which 
downstream flow conditions are given, m. See fig. 7 and sec- 
tion (e), p. 27. 


RHOUT, RTOUT r-Coordinates corresponding to ZHOUT and ZTOUT, m. Leave 

these spaces blank if LROT = 0. If RADOUT is given as input 
(LSFR = 1) , RHOUT and RTOUT cannot be equal to each other. 


SFOUT 


RADOUT 


PROP 


LOSOUT 


Array of values of stream function for input points from hub to 
shroud along line on which downstream flow conditions are given. 
SFOUT is given when LSFR = 0. See fig. 7 and section (e), p. 27. 

Array of r-coordinates of input points along line from hub to 
shroud on which downstream flow conditions are given, m. 
RADOUT is given when LSFR=1. Seefig. 7 and section (e), p. 27. 

Array of absolute total pressures p^ at input points along line 
from hub to shroud on which downstream flow conditions are 

o 

given, N/m . PROP is given when LTPL = 0. See fig. 7 and 
section (e), p. 27. 

Array of fraction of absolute total pressure loss (p^ - p^)/ 

Pq at input points along line from hub to shroud on which 
downstream flow conditions are given. LOSOUT is given when 
LTPL = 1. See fig. 7 and section (e), p. 27. 
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For the case without blades, omit input cards for the variables LAMOUT through 

BETATE, inclusive. 

LAMOUT Array of values of absolute whirl ^rV^ at i n P u t points along line 

from hub to shroud on which downstream flow conditions are 
given, m 2 /sec. LAMOUT is given when LAMVT = 0. See fig. 7 
and section (e), p. 27. 

VTHOUT Array of values of absolute tangential velocity ( VA at input 

' 'o 

points along line from hub to shroud on which downstream flow 
conditions are given, m/ sec. VTHOUT is given when 
LAMVT = 1. See fig. 7 and section (e), p. 27. 

ZBL Two-dimensional array of z-coordinates of points describing 

mean blade surface, m. See figs. 6 and 10 and section (f), 
p. 28. This surface is described by a series (from 2 to 50) of 
blade sections from hub to shroud. The hub section is given 
first, followed by successive sections up to the shroud. 

RBL Two-dimensional array of r- coordinates, corresponding to ZBL, 

of points describing mean blade surface, m. See fig. 6 and sec- 
tion (f) , p. 28 . 

THBL Two-dimensional array of 9 - coordinates, corresponding to ZBL, 

of points describing mean blade surface, rad. See fig. 10 and 
section (f) , p. 28. The tangential coordinate 9 is positive in 
direction of positive rotation (fig. 4). The origin of 
0-coordinates can be anywhere around the circumference. 

THBL is given only when LB LAD is 0 or 1. 

TNBL Two-dimensional array of blade thicknesses normal to mean cam- 

ber line, corresponding to ZBL, RBL coordinates, m. When 
there is little blade lean, small blade surface curvatures, and 
near-parallel suction and pressure blade surfaces, it is feasible 
to use the blade normal thicknesses for TNBL array (LBLAD = 0). 
Otherwise it is recommended to give either the tangential thick- 
ness, TTBL (LBLAD = 1), or the suction and pressure blade sur- 
face coordinates, TH1BL and TH2BL (LBLAD = 2). See fig. 10(a) 
and section (f), p. 28. 

TTBL Two-dimensional array of blade tangential (circumferential) 

thicknesses, corresponding to ZBL, RBL coordinates, rad. 

TTBL is the blade tangential thickness in meters, divided by 
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TH1BL 

TH2BL 

BETALE 

BETATE 

ZHST 

RUST 

ZTST 

RTST 

FLFR 

PERCRD 


RBL. TTBL is given only when LB LAD = 1 . See fig. 10(b) and 
section (f), p. 28. 

Two-dimensional array of upper blade surface coordinates (fig. 
(10(c)) corresponding to ZBL, RBL coordinates, rad. TH1BL is 
given only when LBLAD = 2. See fig. 10(c) and section (f) , p. 28. 

Two-dimensional array of lower blade surface coordinates (fig. 
(10(c)) corresponding to ZBL, RBL coordinates, rad. TH2BL is 
given only when LBLAD = 2. See fig. 10(c) and section (f) , p. 28. 

Array of leading- edge, mean- camber-line tangency angle, deg. 
BETALE is given only when LETEAN = 1; omit this card if 
LETEAN = 0. See fig. 10(d) and section (f), p. 28. 

Array of trailing- edge, mean- camber- line tangency angles, deg. 
BETATE is given only when LETEAN = 1; omit this card if 
LETEAN = 0. See fig. 10(d) and section (f), p. 28. 

Array of z-coordinates of intersections of hub-shroud output sta- 
tion lines with hub profile, m. Omit this card if NOSTAT = 0. 

See fig. 9 and section (h), p. 29. 

Array of r-coordinates corresponding to ZHST array. Omit this 
card if NOSTAT = 0 or if LROT = 0. See fig. 9 and section (h), p. 29. 

Array of z-coordinates of intersections of hub- shroud output sta- 
tion lines with shroud profile, m. Omit this card if NOSTAT = 0. 
See fig. 9 and section (h), p. 29. 

Array of r-coordinates corresponding to ZTST array. Omit this 
card if NOSTAT = Oor if LROT = 0. See fig. 9 and section (h), p. 29. 

Array of values of stream function designating streamlines along 
which output is to be printed. Omit this card if NSL = 0. See 
fig. 9 and section (h), p. 29. If no cards are given for FLFR 
array (i.e. , NSL = 0), the program will automatically assign the 
following 11 values to FLFR: 0.0, 0.1, 0.2, ..., 0.9, 1.0. 

Array of streamwise fractional distances corresponding to the 
PERLOS array. Omit this card if NLOSS = 0. PERCRD should 
have 0. and 1. for the first and last values, respectively. When 
there are blades, PERCRD is the fractional distance between the 
leading and trailing edges. When there are no blades, PERCRD 
is the fractional distance between the upstream and downstream 
input stations. 
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PERLOS 


Array of fractional loss distribution in the streamwise direction 
within the blade, or between upstream and downstream input sta- 
tions when there are no blades. Omit this card if NLOSS = 0. 
PERLOS should have increasing values from 0.1 to 1. , although 
the final value need not reach 1. 


The remaining seven integer variables, beginning with IMESH, are used to indicate 
what output is desired. The program reduces the input mass flow by REDFAC and 
solves the resulting problem iteratively with a stream-function analysis. The results 
of this analysis are then used with full mass flow to obtain a transonic analysis by 
velocity- gradient methods. Output can be obtained after each reduced- mass- flow iter- 
ation and also after the final transonic solution. 

For all these variables except ISUPER, the integer given indicates the multiple of 
iterations at which the user wishes the output associated with the variables to be 
printed or plotted. A zero in any of these indicates that the output associated with that 
variable is to be omitted. A 1 will print or plot the output on every iteration. A 3, for 
example, would give the output on every third iteration. Any nonzero integer will cause 
output to be given on the first and last iterations and after the transonic solution, in 
addition to the output called for at other iterations. ITSON, however, does not give 
output on the first iteration. A large integer will obviously give only the first iteration 
and the final converged iteration. Care should be used not to call for more output than 
is really useful. The following list gives the output associated with each of these 
variables: 


IMESH 


ISLINE 


ISTATL 


Major output at every mesh point of the orthogonal mesh, that is, mesh 
point indices and coordinates, stream-function value, relative velocity 
components and magnitude, critical velocity ratio, and flow angles. This 
output is usually not requested, except for debugging purposes, since 
output called for by ISLINE or ISTATL is preferable in most cases. 

Major output along streamlines (indicated by values in FLFR) at each point 
where streamlines are crossed by vertical mesh lines of the orthogonal 
mesh. Output includes z, r, and m streamline coordinates, relative 
velocity components and magnitude, critical velocity ratio, flow angles, 
and streamline curvature. Where the streamline passes within the blade 
region, an estimate of suction- and pressure-surface velocities is also 
printed. 

Major output along station lines from hub to shroud, at locations specified 
by ZHST and ZTST arrays and at values of stream function specified by 
FLFR array. Output corresponds to that given for ISLINE with the addi- 
tion that absolute velocity components, as well as density, pressure, and 


24 



IPLOT 


ISUPER 


ITSON 


IDEBUG 


temperature for static conditions and for absolute and relative stagnation 
conditions are printed. ISTATL should be zero if NOSTAT is zero. 

Plotting indicator requesting output to be plotted on microfilm. Any non- 
zero value in IPLOT will cause the input data and generated orthogonal 
mesh to be plotted. Also, at each iteration that is a multiple of IPLOT, 
streamlines will be plotted and meridional and surface velocities will be 
plotted for each streamline value. These will also be plotted after the 
final transonic solution. 

Integer (0 or 1) indicating which solution (subsonic or supersonic) of the 
velocity-gradient equation is desired. If ISUPER = 0, only the subsonic 
solution will be printed. If ISUPE R = 1, both subsonic and supersonic 
solutions will be printed. 

Integer indicating when information is desired for use in calculating input 
for the TSONIC program (ref. 3). If ITSON = 0, no information will be 
given for TSONIC. Otherwise, information will be listed with other out>- 
put. Usually, the user only wants TSONIC data after subsonic conver- 
gence is reached. Using a high value for ITSON (i. e. , >20) will achieve 
this result. 

Integer indicating whether additional debug output is desired. If 
IDEBUG = 0, no extra output is printed. If IDEBUG > 0, blade geometry 
information on the alternate mesh and the coefficients of the finite- 
difference equations, as well as 21 arrays of debug output on the or- 
thogonal mesh, are printed. Fifteen of these arrays change and are re- 
printed after each iteration of the reduced-mass-flow solution that is a 
multiple of IDEBUG. 


Special Instructions for Preparing Input 

It is unusual to have no errors in input to MERIDL the first time any new data set is 
run. Therefore, input should be checked thoroughly before it is submitted. Errors are 
commonly made for the following reasons: inconsistent units; improper sign on co, V^, 
or whirl, input for arrays not agreeing with the input bounds for those arrays; and up- 
stream and downstream input not being of the form specified by LSFR, LTPL, and 
LAMVT. Also geometric input into the hub and shroud arrays and the blade- geometry 
arrays should be smooth enough that the hub, shroud, and blade sections will be fit well 
with cubic spline curves (see section (g), p. 29). These input geometric arrays are all 
plotted on microfilm by the program after spline fitting is completed; the microfilm 
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output will indicate whether or not the input was smooth. All output should be checked, 
especially that from a new input data set, to see if it is reasonable. 

(a) Units of measurement. - The International System of Units (ref. 14) is used 
throughout this report. However, the program does not use any constants that depend 
on the system of units being used. Therefore, any consistent set of units may be used 
in preparing input for the program. For example, if force, length, temperature, and 
time are chosen independently, mass units are obtained from Force = Mass x Acceler- 
ation. The gas constant R must then have the units of (Force x Length)/ (Mass x 
Temperature). Density is mass per unit volume, and mass flow is mass per unit time. 
Output then gives velocity in the chosen units of length per unit time. Since any consis- 
tent set of units can be employed, the output is not labeled with any units. 

(b) Damping factors FNEW and DNEW . - The input variables FNEW and DNEW are 
used as damping factors on the F fc and 9(rV 0 )/9t terms of equation (Al) . During each 
iteration of the subsonic or finite- difference portion of the program, new values for F^. 
and 9(rVg)/9t are calculated. The calculated changes to these terms are often so 
large that, if the full change was accepted on each iteration, the solution would diverge. 
Using values of FNEW and DNEW less than 1. 0 allows only a portion of these changes to 
be used. The value of FNEW or DNEW is the fraction of the predicted change in F^. or 
9(rV 0 )/9t to be added to the previous values of these variables. Therefore, if F^. or 
9(rV 0 )/9t does not converge, FNEW or DNEW should be reduced. However, this will 
reduce the rate of convergence so that caution should be used that FNEW and DNEW are 
not made too small. Because of the reduced rate of convergence, VELTOL is auto- 
matically reduced to obtain the same accuracy as when FNEW and DNEW are 1.0. 

Limited experience has shown that for subsonic flows, where the gradients of flow 
variables are not large, converged solutions can be obtained with FNEW and DNEW 
equal to 1. 0. For most practical flows, however, this is not the case and some damp- 
ing is required. Values of FNEW and DNEW in the range of 0. 5 are most commonly 
used, although at times values as low as 0. 05 have been used. (Note that FNEW and 
DNEW do not have to be equal to each other.) 

The user will have to gain experience with FNEW and DNEW before he can use 
them effectively to maximize rates of convergence. Maximum and minimum values and 
the maximum change for F^ and 9(rV 0 )/9t (DVTHDR) are printed with the output on 
each iteration. The maximum relative change in velocity at any point is also printed; 
it should gradually approach VELTOL for convergence. By observing these values, the 
user can learn the effects of changes in FNEW and DNEW for different runs of the pro- 
gram. 

(c) Hub and shroud flow- channel geometry . - The hub and shroud flow-channel 
geometry is specified in the ZHUB, RHUB and ZTIP, RTIP arrays. Both of these curves 
must have the same z-origin (typically the blade leading edge at the hub) . These two 
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arrays must extend far enough upstream and downstream to cover the upstream and 
downstream boundaries of the orthogonal mesh, as well as the upstream and down- 
stream input stations where streamflow data are given. If they do not extend this far, 
they will be linearly extrapolated and an incorrect flow channel may result. Relatively 
few points are needed to describe these smooth surfaces (2 to 10 is a typical range for 
NHUB and NTIP) in order to have the program calculate smooth, accurate spline fits of 
these surfaces (fig. 6). 

If the user knows the amount of the boundary- layer blockage along the hub and 
shroud profiles, he should revise the hub and shroud arrays to include this blockage 
effect. In this situation the output will represent a more realistic flow condition. 

(d) Orthogonal mesh . - The number of orthogonal mesh lines is specified by MBI, 
MBO, MM, and MHT; and the positioning of the mesh is specified by ZOMIN, ZOMBI, 
ZOMBO, and ZOMOUT. When ANGROT is specified (LROT = 1), ROMIN, ROMBI, 
ROMBO, and ROMOUT must also be specified. The four z- coordinates, all located on 
the hub, must use the same z-origin as all other geometric input. The ZOMBI and 
ZOMBO locations are usually close to the leading and trailing edges of the blade at the 
hub, although they do not have to correspond to these locations (fig. 8). Mesh-size 
spacing in the horizontal direction is established by the relation between the four 
z-coordinates (ZOMIN, etc.) and the number of mesh spaces requested. For instance, 
MBO - MBI evenly spaced mesh spaces along the hub will be located in the distance 
given by ZOMBO - ZOMBI. Mesh spacing in the vertical direction is determined by 
locating MHT spaces in the hub- shroud distance. Usually, none of the vertical mesh 
lines follow the blade leading- or trailing-edge lines since the blade edges are not 
usually orthogonal to the hub or the shroud (fig. 8) . So some vertical mesh lines will 
pass in or out of the blade region as they progress from hub to shroud. The suggested 
number of mesh lines between ZOMBI and ZOMBO that are required to cover the blade 
adequately is 15 to 30, depending on blade geometry. 

(e) Upstream and downstream flow conditions . - Upstream flow conditions are 
given in the SFIN (or RADIN), TIP, PRIP, and LA MIN (or VTHIN) arrays, which are all 
of length NIN. Downstream conditions are given in the SFOUT (or RADOUT), PROP 

(or LOSOUT), and LAMOUT (or VTHOUT) arrays, which are all of length NOUT 
(fig. 7). Upstream and downstream flow conditions are used, along with the assump- 
tion of conservation of angular momentum along streamlines, to establish the flow field 
in the upstream and downstream portions of the solution region. 

Ordinarily, the upstream and downstream flow conditions are given as a function of 
stream function (SFIN and SFOUT). In this case these input values apply at all points 
along streamlines upstream or downstream of the blade. If, in addition, the whirl is 
specified (LAMVT = 0) and there are blades (MBI ^ 0), the ZHIN, ZTIN and 
ZHOUT, ZTOUT inputs are superfluous, so that these variables need not be specified. 
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However, in most cases, legitimate values must be supplied for ZHIN, ZTIN and 
ZHOUT, ZTOUT. In this case the upstream conditions are given on a straight line that 
passes through the two points given by ZHIN on the hub and ZTIN on the shroud. Down- 
stream conditions are given on a straight line that passes through ZHOUT on the hub 
and ZTOUT on the shroud. These lines may lie anywhere in the regions from the blade 
edges upstream and downstream to the boundaries of the orthogonal mesh (figs. 7 and 8). 
If LROT is 1, values for RHIN, RTIN, RHOUT, and RTOUT must also be given as input. 

The arrays of upstream and downstream input do not necessarily have to extend all 
the way from the hub to the shroud or lie on radial lines. They will be linearly extrap- 
olated to the hub and the shroud, if necessary, by the program, should the user only 
give data in a portion of the flow channel. 

(f) Mean blade surface and thickness coordinates. - The blade shape is described 
from hub to shroud by four arrays, all of which are two dimensional (figs. 6 and 10). 
Each of these arrays has NBLPL blade sections or planes, with NPPP points in each of 
these sections. When giving data for each of these four arrays, start each new section 
of data (NPPP points) at the beginning of a new card. There does not have to be any 
geometric relation between analogous points on adjacent blade sections. All the ZBL 
data for all the sections are given, followed by all the RBL data, etc. The origin for 
the z-coordinates of ZBL should be the same as that used for all other z-coordinate in- 
put arrays. Of the four arrays necessary to describe the blade shape, ZBL and RBL 
are always given. There are three options, controlled by LBLAD, for the other two 
arrays (see LBLAD definition in the input dictionary). 

When LBLAD is 0 or 1, the THBL array must be given (figs. 10(a) and (b)). The 
THBL array is for input 0 - coordinates of the blade mean camber surface. Alterna- 
tively, the midchannel stream surface, obtained from a blade-to-blade solution (ref. 3), 
may be used for the THBL array. The program provides a solution on the midchannel 
stream surface. It obtains that stream surface by smoothly fairing the inlet and out- 
let flow angles onto the input blade surface at an appropriate distance back on the blades 
from the leading and trailing edges. This distance is a function of solidity (appendix F) . 

The TNBL array is given when LBLAD = 0 (fig. 10(a)). The TNBL array is for in- 
put blade thicknesses normal to the blade-section mean camber line and lies on a sur- 
face of revolution cutting through the blade (figs. 6 and 10). So, in general, the thick- 
nesses lie on a curved line whose ends may be at different radii and may or may not be 
normal to the blade surfaces (fig. 10(a)). Because it is difficult to obtain proper block- 
age or proper TSONIC input, caution should be exercised when giving TNBL input. Al- 
so, when TNBL is given as input, the THBL coordinates must be for the blade mean 
camber surface and not the midchannel stream surface. 

The TTBL array is used for input tangential (circumferential) blade thicknesses 
(when LBLAD = 1) and is given in radians (i.e. , t^/r) (fig. 10(b)). Thick blades, 
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blades with high curvature, or blades with significant lean should use TTBL input in 
preference to TNBL input. 

The TH1BL and TH2BL arrays are used to give blade surface 0-coordinates (when 
LB LAD = 2) (fig. 10(c)) . This input is similar to that used in the TSONIC program 
(ref. 3). 

To obtain a smooth solution near the leading and trailing edges, it is recommended 
to use nonzero values of leading- edge and trailing- edge thicknesses, as figure 10 indi- 
cates. 

In addition to the four arrays necessary to describe the blade, there are two one- 
dimensional arrays (BETALE and BETATE) that are optional. These arrays specify 
the leading-edge and trailing-edge mean-camber-line tangency angle in degrees. The 
angles are measured on the input blade-section surfaces (fig. 10(d)). Use LETEAN = 1 
when BETALE and BETATE are given, and LETEAN = 0 when they are not given. 

The first blade section given at the hub or the last one at the blade tip does not 
necessarily have to conform to the hub or shroud profile. It can be given within the 
flow region, crossing the boundary, or completely outside the boundary (fig. 6). Ex- 
trapolation or interpolation will be used when necessary to obtain blade data where the 
blades meet the hub and shroud profiles. 

The user should attempt to give smooth data for these four blade- shape arrays. 

The program makes use of spline fits of these curves to compute gradients that are 
used throughout the program. The mean- camber-line first and second derivatives, as 
well as microfilm plots of the blade sections, given as output will indicate whether input 
data for these arrays were acceptable. 

(g) How to specify points for spline curves . - All the input arrays are fit with cubic 
spline curves for the purposes of interpolating, calculating derivatives, or any other 
required calculation. A cubic spline curve is a piecewise cubic polynomial that ex- 
presses mathematically the shape taken by an idealized spline passing through the given 
points. Reference 15 describes the method used for determining the equation of the 
spline curve. By this method, smooth curves can be specified accurately with a few 
points, usually not more than four or five. Curves with uneven places, dips, or highly 
variable curvatures require more points and are more difficult to fit properly. As a 
guide, enough points should be specified so that a physical spline passing through these 
points would accurately follow the curve. The minimum number of points to follow the 
curve should be used, since the closer the spline points are, the greater the effect of an 
inaccuracy in a coordinate. 

(h) Requests for output data . - The seven variables IMESH, ISLINE, ISTATL, 

IP LOT, ISUPER, ITSON, and IDEBUG all request different portions of output from the 
program. The optional arrays ZHST, RHST, ZTST, RTST, and FLFR are also used in 
this regard. The user should be careful to request only the output needed. Usually, 
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Shroud 



cross user-designated station lines 

Streamlines and station lines 

Orthogonal mesh lines 

Figure 11. - Location of three major types of output. 


only one of the three major types of output (at the mesh points (IMESH), along stream- 
lines (ISLINE), or along hub- shroud stations (ISTATL)) is needed, since the ISLINE and 
ISTATL outputs are obtained from calculated mesh-point data by interpolation (fig. 11). 
Most likely, streamline output (ISLINE) or station-line output (ISTATL) would be chosen. 
The frequency at which output is requested is also significant; this is controlled by the 
values given to IMESH, ISLINE, and ISTATL. 

ISTATL output should be requested only if NOSTAT is greater than zero (and the 
ZHST, RHST, ZTST, and RTST arrays are given). If LROT = 0, omit the RHST and 
RTST cards. Output will be given along straight lines from hub to shroud that connect 
corresponding points in these arrays. Through the values in these arrays, the user 
can control exactly the locations at which his output is given. These lines do not have 
to be radial. Output can be requested at several upstream and downstream locations, 
at the blade leading and trailing edges, and at several stations within the blade at per- 
centages of chord from the leading to the trailing edge. To obtain output at the leading 
and trailing edges, the values in ZHST, RHST, ZTST, and RTST should correspond 
(within a tolerance of 1 percent of local chord) to the intersection points of the leading 
and trailing edges with the hub and shroud. 

The array FLFR need be given only if the user wishes to specify the values of 
streamline along which he wishes output (when ISLINE or ISTATL are used) . If the 
user does not specify otherwise (using FLFR), the program will automatically give out- 
put at 10-percent streamlines. 
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(i) Incompressible flow . - Incompressible flow cases can be executed on MERIDL 
by using special input for GAM and AR. Use GAM = 0. and set AR equal to the fluid 
density (RHO). In this case, TIP, PRIP, and PROP are not used by the program, but 
blank cards should be furnished with the input. No other special input is required for 
an incompressible case, except that REDFAC should be set to 1.0 to avoid the transonic 
solution, which would be redundant and less accurate. 

(j) Straight infinite cascade . - The program is primarily designed for circular, 
stationary, or rotating blade rows; but the input can be adapted to apply to a straight 
infinite cascade as well. Since the radius for such a cascade would be infinite, an arti- 
ficial convention must be adopted. The user should pick a large average radius for his 
cascade, in the neighborhood of r= 1000. Then, since the blade pitch P is known, the 
approximate number of blades to use can be calculated from NBL = (27rr)/P. This value 
must be rounded to an integer for use in the input. With this value of NBL and a large 
mean radius, all the r- coordinate arrays of input (RHUB, RTIP, RADIN, RADOUT, and 
RBL) can be established about the mean radius, and a value of MSFL can be calculated. 
Using this artificial input at large radius, the program will obtain a solution for a pitch 
that varies very slightly from hub to shroud and, therefore, simulates almost exactly a 
straight infinite cascade. 

(k) Choosing a value for REDFAC . - If possible, REDFAC should be 1.0. However, 
if there is locally supersonic flow, REDFAC must be less than 1. 0, usually between 0. 5 
and 0. 95. For the best accuracy the largest possible value of REDFAC should be used 
(see the section NUMERICAL EXAMPLE, p. 70). When REDFAC is 1.0, the full mass 
flow will be used in the finite-difference solution for the stream function, and no tran- 
sonic velocity- gradient calculation will be made. 

( l ) Choosing a value for ANGROT . - If the flow is close to axial (within about 45° 
from axial), ANGROT need not be given (use LROT = 0). If the flow deviates more than 
45° from axial, use LROT = 1 and specify ANGROT in degrees so as to minimize the 
maximum slope of the hub or shroud from the rotated axis. For example, an impeller 
with axial inlet and radial discharge should have ANGROT = 45 °. A radial inflow tur- 
bine rotor with axial discharge should have ANGROT = - 45 °. 

OUTPUT 

There are four different types of output generated by the MERIDL program: 

(1) Main output - controlled by the variables IMESH, ISLINE, ISTATL, and ISUPER 

(2) Debug output - controlled by IDEBUG 

(3) Information for the TSO NIC program - controlled by ITSON 

(4) Plotted output - controlled by IPLOT 
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Most of this output is optional and is controlled by the final input card, as already de- 
scribed. 

The output controlled by each of the variables IMESH, ISLINE, and ISTATL is es- 
sentially the same but is given at different locations for the convenience of the user. 

The IMESH output is given at the orthogonal mesh points along horizontal mesh lines, 
as indicated in figure 11. The ISLINE output is given along streamlines where the 
streamlines are intersected by the vertical orthogonal mesh lines (fig. 11). The ISTATL 
output is given from hub to shroud along station lines (fig. 11) where these lines are 
intersected by the streamlines. Also, additional ISTATL output is given beyond that 
given with IMESH or ISLINE output. 

In the following sections, output is presented from the problem solved in the section 
NUMERICAL EXAMPLE (p. 70). Since the complete output would be lengthy, only the 
first few lines of each section of output are reproduced here. For debug output 
(IDEBUG ^ 0), output labels are simply internal variable names, which are defined in 
the main dictionary of part II (ref. 6) . 

The following three sections discuss the different sections of printed output, the 
plotted output, and all possible error messages. 


Printed Output 

Table I presents the printed output from the numerical example. Each section of 
this output has been numbered to correspond to the following description: 

(1) The first output is a listing of the input data. Variable names are used as 
labels, and the output corresponds to the input form (fig. 5). This output is listed for 
every run, regardless of the values given to IMESH, etc. 

(2) The second output gives information to assist the user in checking input blade 
geometry. This output is listed for every run. For each input blade section from hub 
to shroud the following is given: ZBL, RBL, THBL, and TTBL, all defined as in the 
input dictionary; blockage at each input point (tangential blade thickness TTBL divided 
by blade- to- blade pitch in radians); SZRBL, the meridional arc length of the input blade 
section; and DTHDSP and DTHDSP2, first and second derivatives of THBL against 
SZRBL along the input blade sections. The smoothness of DTHDSP and DTHDSP2 indi- 
cates where inaccurate data or errors may have occurred in input blade geometry. 

(3) The third output corresponds to IDEBUG. It has four principal sections, as the 
printed output indicates. The first section gives blade geometry where the blade is in- 
tersected by the 21x21 alternate mesh used for calculating gradients of the mean 
camber surface. The z, r, 9, and tangential thickness coordinates are given, followed 
by the gradients of 0, calculated by spline curves, along the two directions of the alter- 
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TABLE I. - SAMPLE OUTPUT 
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STREAM FUNCTION, INTERIOR VELOCITIES, VELOCITY COMPONENTS, AND ANGLES 
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STREAM PUNCTICN, INTERIOR VELOCITIES, VELOCITY COMPONENTS, ANGLES, AND SURFACE VELOCITIES 
ALONG LINES FROM HUB TO SHROUD AT VARIOUS STATIONS THROUGH THE BLADE ROM 
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nate mesh. This first section of debug output is given only once. The second section 
of debug output lists orthogonal mesh quantities that are held constant. These quantities 
are also listed only once. The third and fourth sections of debug output list the coeffi- 
cients of the matrix equation that is solved for the stream function and the orthogonal 
mesh quantities that change with each iteration. These two sections are listed follow- 
ing an iteration of the solution, as called for by the value of IDEBUG. The output 
giving the changing orthogonal mesh quantities contains two different sets of variables. 

(4) The fourth output indicates the calculated value of overrelaxation factor ORF to 
be used in the successive overrelaxation solution of the finite- difference equations. 

This matrix equation for the stream function u is solved iteratively (this is referred to 
as the "inner iteration") during each of the "outer iterations" of the reduced-mass-flow 
solution. The ORF output is given automatically on each run of the program. 

(5) The fifth output lists the calculated LOSOUT array at the points along the line 
where input downstream flow conditions are given. This array is printed if PROP is 
given as input. 

(6) The sixth output lists the ILE and ITE arrays. For the J th horizontal mesh line 
from the hub, ILE(J) and ITE(J) are the vertical mesh line numbers of the mesh points 
that are within the blades and closest to the leading and trailing edges. 

(7) The seventh output gives convergence information about DVTHDT, FT, and the 
relative velocity. Convergence is determined by the maximum relative change in veloc- 
ity. When the maximum change becomes less than the input value of VELTOL (reduced 
by the minimum of FNEW or DNEW) , the stream-function (either reduced or full mass 
flow) solution is considered converged. Also given are the indices where the maximum 
change occurred, the average relative change in velocity, and the number of uncon- 
verged points . 

If there is a convergence problem, the behavior of DVTHDT and FT are of interest. 
Therefore, the maximum and minimum values and the maximum change in DVTHDT 
(9(rV^)/3t) and FT(F^) are printed for each iteration. The input variables DNEW and 
FNEW are used to control the fraction of the calculated changes of these variables to be 
used for the next following iteration. If the maximum change in either DVTHDT or FT 
is increasing with each iteration, it is suggested that either DNEW or FNEW be de- 
creased. 

(8) The eighth output corresponds to IMESH for one of the iterations of the stream- 
function solution. The output is given along each horizontal mesh line at each mesh 
point. The output given at each mesh point includes the following: z- and r- coordinates 
of the point, stream function u, meridional velocity W , relative tangential velocity 
Wg, relative velocity W, critical velocity ratio W/W cr , meridional flow angle a, 
relative flow angle (3, and mesh angle cp. 


48 



(9) The ninth output corresponds to ISLINE for one of the iterations of the stream- 
function solution. This output is given along each streamline, corresponding to a given 
stream- function value. The points along the streamline correspond to where it is inter- 
sected by the vertical mesh lines. The origin for the m-coordinate is located where 
each streamline begins at the upstream boundary of the orthogonal mesh. The output 
given at each streamline point includes the following: z-, r-, and m- coordinates of the 
point, meridional velocity relative tangential velocity W^; relative velocity W; 

critical velocity ratio W/W^; meridional flow angle a; relative flow angle /3; 
streamline curvature 1/ r c ; and within the blade- to- blade passage, the estimated 
suction-surface and pressure-surface velocities. 

(10) The 10th output corresponds to ISTATL for one of the iterations of the reduced- 
mass-flow solution. Instead of output being given along horizontal mesh lines (IMESH) 
or streamlines (ISLINE) , it is now given in the other direction, along lines from hub to 
shroud. This output is given along each of the station lines specified by the input 
ZHST, ZTST arrays. Output is given at each point where these station lines are crossed 
by the streamlines (the values in the FLFR array) . The output given at each point is 
identical to the output given for ISLINE, with the addition of the stream function u. The 
second row of output includes the following additional variables: a repeat of z- and 



Figure 12. - Definition of incidence angles. 
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r- coordinates of the point, static pressure p, static temperature T, static density p, 
absolute tangential velocity V^, absolute velocity V, absolute total pressure p\ abso- 
lute total temperature T', absolute flow angle /3 ^ g , relative total pressure p", and 
relative total temperature T". 

(11) The 11th output consists of blocked and unblocked incidence and deviation angles 
and blade angles at the blade leading and trailing edges where each of the horizontal 
mesh lines intersects the blade. The blocked incidence i^ or deviation is based on 
the velocity diagram within the blade (subscript bf in fig. 12); the unblocked inci- 
dence i| g is based on the free-stream velocity diagrams (subscript fs in fig. 12) . 

This output is printed after each iteration in which output is given for either IMESH, 
ISLINE, or ISTATL. These angles are defined as shown in figure 12. 

(12) The 12 th output corresponds to ITSON. It is a listing of the information re- 
quired to prepare input for the TSONIC blade- to-blade analysis program of reference 3. 
(The TSONIC program is being updated and will be reported in the future. The defini- 
tions of all the input variables for TSONIC will be given in the updated version of 
TSONIC.) The TSONIC program was run for both the hub and tip blade sections of the 
compressor blade used in the numerical example of this report. The input numbers 
shown in the input form in figure 13 are those computed for use with the TSONIC run at 
the hub. 


The information printed by MERIDL can in some cases serve as input for TSONIC 
without revision. In other cases, some of the input may need revision. The informa- 
tion from spline fits (derivative, second derivative, and curvature) of the THSP- 
against-ZMSP arrays for each blade surface is printed alongside these arrays. The 
user may wish to make small changes in the BETI and BETO tangency angles, which 
serve as end conditions for these spline fits, in order to smooth the resultant first and 
second derivatives. 


The final five arrays of stream- channel data - ZMRSP, RMSP, BESP, WOWCR, 
and PLOSS - can in general be used as they are printed, unless the user wishes to re- 
duce the number of points (NRSP) in them. The BESP array is calculated to correspond 
to the given mass flow ZMSFL. Any value could be used for ZMSFL if BESP is chosen 
accordingly. The value of ZMSFL calculated by MERIDL is 1 percent of the mass flow 
between any two blades. The BESP array is sometimes bumpy and may also have to be 
smoothed before it can be used as input to TSONIC. The derivatives of BESP against 
ZMRSP are printed and will indicate if smoothing is necessary. The WOWCR array is 
only needed as TSONIC input if REDFAC < 1.0 is used for TSONIC. The printed PLOSS 
array (NLOSS > 1) can be used for TSONIC rather than a single value for PLOSS if a 
nonlinear loss distribution through the blade is used in MERIDL. 

(13) The 13th output is analogous to outputs 8, 9, 10, 11, and 12 for the variables 
IMESH, ISLINE, ISTATL, and ITSON but is given after the transonic velocity-gradient 
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Figure 13. -TSONIC input form. (Data shown are for the numerical example at the hub. ) 
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solution. Two solutions can be obtained by the velocity- gradient method: the larger or 
"supersonic" and the smaller or "subsonic." The input variable ISUPER controls which 
of these is obtained. The solution listed here is the smaller (i.e., the subsonic solu- 
tion). If REDFAC = 1.0, no velocity- gradient solution will be obtained or printed be- 
cause the final iteration of the stream-function solution is the correct solution to the 
problem. 


Plotted Output 


Since the printed output from a MERIDL run can be quite large, microfilm plots 
are made to enable the user to quickly check the quality of both his input and output 
data. The coding that generates these plots uses NASA Lewis in-house plot routines 
and would have to be recoded for operation on another system. The principal sections 
of input data plotted are the upstream and downstream distributions of flow properties; 
the input blade sections; and the hub, shroud, and blade leading- and trailing- edge 
geometry. The frequency at which output data are plotted is controlled by the variable 
IPLOT. The plotted output data begin with the generated orthogonal mesh. Then for 
each iteration of the solution (indicated by IPLOT) and for the transonic solutions, 
streamlines are plotted, as well as midchannel and blade surface velocities along each 
streamline from hub to shroud. The user should carefully check the plots of his input 
data to ensure that the program's spline fits of these data are smooth. 

Selected examples of some of the microfilm plots generated by the numerical ex- 
ample are presented in figure 14. The description of these plots is as follows: 

(1) Figures 14(a) to (e) present the input upstream and downstream flow conditions. 


These quantities ^T!, pj, X or (Vg) , p^ or total pressure loss, (rV^ or^V^ 

are plotted against either stream function or radius, whichever was used as input. 

(2) Figure 14(f) shows the optional input of total-pressure-loss distribution against 
fraction of chord. 

(3) Figures 14(g) and (h) indicate hub and shroud blade sections plotted from input 


data supplied by the user. All blade sections given as input are plotted. 

(4) Figure 14 (i) shows the hub and shroud channel geometry and the blade leading 
and trailing edges. 

(5) Figure 14 (j) shows the generated orthogonal mesh. 

(6) Figure I4(k) shows the streamline pattern for the final iteration of the reduced- 
mass-flow solution. 
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DIRECTION 



Z DIRECTION 

(j) Orthogonal mesh in meridional plane. 
Figure 14. - Continued. 
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(7) Figures 14 (l) and (m) show meridional velocities and blade surface velocities 
for the hub and shroud blade sections for the same iteration. A similar plot is made 
along each of the streamlines from hub to shroud. 

Streamline and velocity plots similar to figures 14(k), ( l ), and (m) are made after 
each iteration (indicated by IPLOT) of the finite-difference, stream-function solution 
and also after the velocity- gradient solution. 


Error Messages 

A number of error messages have been incorporated into the program. These 
error messages are listed here, followed by suggestions for finding and correcting the 
cause of the error: 

MM, MHT , NHUB, NT IP , N I N , N OUT , N B LPL , NPPP, NOSTAT, NSL.NLOSS, LSFK , LTPL , LA PI VT , LROT , LBLA D , OR LETEAN 
IS TCC LARGE OR TOO SHALL 

The input dictionary gives the maximum and minimum values for all these input varia- 
bles. 


WHEN UPSTREAM AND DOWNSTREAM INPUT DATA ARE GIVEN AS A FUNCTION OF RADIUS, 

THERE MUST BE A CHANGE IN VALUE BETWEEN RHIN AND RTIN AND ALSO BETWEEN PHOUT AND RTOUT 
AND A CORRESPONDING CHANGE IN THE RADIN AND RADOUT ARRAYS 

RHEN must differ from RTIN (and RJHOUT from RTOUT) so that there will be a variation 
in the RADIN (and RADOUT) arrays, which are used to identify points along the hub- 
shroud lines at the upstream and downstream input stations. 

LININT CANNOT FIND INTERPOLATED VALUE 
I = 1 J = 2 

This message should occur only with erroneous geometry input. 

INPUT PROP VALUES ARE LARGER THAN IDEAL TOTAL PRESSURE, RESULTING IN NEGATIVE LOSS 

This message is for information only; the program will continue normally. 

PROGRAM STOPPED IN NEWRHO DUE TO EXCESSIVE STREAM FUNCTION GRADIENT 

This is usually caused by having supersonic flow in some region of the reduced-mass- 
flow solution. Try a smaller value for REDFAC. 
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THE UPSTREAM INPUT WHIRL OR TANGENTIAL VELOCITY IS TOO LARGE 

This message is printed by subroutine NEWEHO or TVELCY if the upstream whirl gives 
a tangential velocity so large as to result in a negative relative stagnation temperature. 
The upstream whirl or tangential velocity given as input should be corrected. 


BEWARE. THERE ARE 1 VERTICAL ORTHOGONAL HE SH LINES WHICH ARE CHOKED. 

LOCATIONS OP THESE LINES APE GIVEN ABOVE AT THE BEGINNING OF THE TRANSONIC OUTPUT. 
OUTPUT ON ANY STATION LINES LOCATED NEAR THESE CHOKED ORTHOGONAL LINES WAY BE IN ERROR 


When this message is printed, the user should be extremely cautious in using station 
line output after the transonic velocity solution. However, there may be little problem 
if the choking mass flow is close to the full mass flow or if the station line is not near 
one of the choked orthogonals. 

A VELOCITY GRADIENT SOLUTION CANNOT BE OBTAINED FOR VERTICAL ORTHOGONAL HESH LINE I - 1 

ANY SUBSEQUENT OUTPUT FOR THAT MESH LINE WAY BE IN ERROR 

If the approximate transonic solution for any particular vertical mesh line cannot be 
found for some reason, this message will be printed. Some of the following error mes- 
sages will be printed to help determine the problem. The program will proceed to the 
next vertical mesh line. 


CONTIN COULD NOT FIND A SOLUTION IN 100 ITERATIONS. 

This message is printed when the transonic flow solution for some vertical mesh line 
cannot be obtained by subroutine CONTIN in 100 iterations. Possible causes for this 
error message could be unusual or impossible upstream or downstream input quantities 
(TIP, PRIP, LAMIN, VTHIN, PROP, LOSOUT, LAMOUT, or VTHOUT arrays). 


ITERATION PROCEDURE HAD TO BE RESTARTED TO AVOID NEGATIVE TEHPEBATURE OR VELOCITY 

HAGNITODE LESS THAN TANGENTIAL VELOCITY, OR AFTER ADJUSTMENT OF STAGNATION TEMPERATURE, 

RESTART OF ITER AT ION * PROC EDUF E (LOOP TO STATEMENT 90) WAS ABORTED AFTER 1000 OR MORE 
TOTAL ITERATIONS (LOOP TO STATEMENT 150). 


This message is printed when the transonic velocity- gradient solution cannot be found 
by subroutine TVELCY in 1000 iterations for some vertical mesh line. Possible causes 
are high input whirl, impossible upstream or downstream flow condition, and complex 
or erroneous blade or wall geometry. 
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THE HAXIHUH HASSFLOW FOB WHICH A SOLUTION COULD BE OBTAINED WAS 0 . 1 2 00 0 C 0E-02 

THE BINIHUH BASSFLOW FOR WHICH A SOLUTION COULD BE OBTAINED WAS 0 . 1 0 1 0000 E-0 2 

THE HAXIHUH VALUE OF W AT THE HUB FOR WHICH A SOLUTION COULD BE OBTAINED WAS 101.1000 

THE BINIHUH VALUE OF W AT THE HUB FOR WHICH A SOLUTION COULD BE OBTAINED WAS 99.79999 


This information is printed whenever the transonic, velocity- gradient solution cannot 
be found for some vertical mesh line. This information has debug value if there is no 
apparent reason for difficulty. 


THE ITERATION PROCEDURE WAS RESTARTED 380 TIHES 
WHUB WAS INCREASED 112 TIHES 
WHUB WAS DECREASED 56 TIHES 

BOUNDARY VALUES (TIPBDT, RHOIP, LAHBDA) WERE ADJUSTED 34 TIHES 


This information is printed whenever the transonic solution cannot be found after re- 
starting the iteration procedure. This information has debug value if there is no 
apparent reason for difficulty. 

HSFL EXCEEDS CHOKING HASS FLOW FOR VERTICAL ORTHOGONAL HESH LINE I = 1 

CHOKING HASS PLOW = 0.950000E-03 

This message is printed if the choking mass flow calculated by subroutine TVELCY is 
less than the input mass flow MSFL. The following message will be printed after all 
vertical lines have been checked: 


CHOKING HASSFLOW IS 0.96000 OF THE INPUT HASSFLOW 
INPUT HASSFLOW = 0.15000E-02 

BINIHUH CHOKING HASSFLOW = 0.1U400E-02 

If the previous message was printed for any vertical mesh line, this message will be 
printed at the end to give the minimum choking mass flow, which is the choking mass 
flow for the blade row estimated by MERIDL. Use of the CHANEL program (ref. 13) 
is recommended to obtain a more accurate estimate of the choking mass flow. 

SUPERSONIC CORRECTION - LEADING EDGE FOR J = 2 

For the full-mass- flow solution when REDFAC < 1, the blockage correction at the lead- 
ing edge requires calculating the density by an iterative procedure. Since the iterative 
procedure depends on whether the flow is subsonic or supersonic, the program first 
checks the subsonic procedure; if this fails, it is usually because the flow is super- 
sonic. If the supersonic solution is found, the preceding message is printed for infor- 
mation only. 
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SUPERSONIC CORRECTION - TRAILING EDGE FOR 3 = 2 


This is the same as the previous message but at the trailing edge. 

NO DENSITY CORRECTION - LEADING EDGE FOR J = 2 

If after the subsonic procedure fails the supersonic procedure is tried and also fails, 
that procedure is aborted and the density within the blade is taken to be the free- stream 
density. This problem would probably occur because of numerical problems associated 
with near-sonic flows, high blockage, and complex or erroneous wall or blade geometry. 
The program will continue after this message, although output near the leading edge 
may be poor. 

NO DENSITY CORRECTION - TRAILING EDGE FOR J * 2 

This is the same as the previous message but at the trailing edge. 


INRSCT HAS FAILED TO CONVERGE IN 20 ITERATIONS 
TOLERANCE - 0.1C1000E-05 

DISTANCE BETWEEN LAST TWO INTERSECTION POINTS = 0.678900E-05 


Subroutine INRSCT finds the intersection coordinates of input station lines, mesh lines, 
and streamlines with the blade leading or trailing edge by an iterative method. If the 
tolerance cannot be met after 20 iterations, the message is printed. If the distance be- 
tween the last two intersection points is only slightly larger than the tolerance, a satis- 
factory solution will be obtained, with some loss of accuracy. If the distance is exces- 
sive, there is probably some error in the geometry input. 


SPLINE ERROR — ONE OF THREE POSSIbLE CAUSES 

1. ADJACENT X POINTS ARE DUPLICATES OF EACH OTHER. 

2. SOME X POINTS ARE OUT OF SEQUENCE. 

3. SOME X POINTS ARE UNDEFINED. 


NUMBER OF POINTS = 3 


X ARRAY 

0.00000 

0.00000 

1.0000 


Y ARRAY 

10.000 

12.000 

15.000 


The spline points for a spline curve must be distinct and given in sequence. If not, the 
spline fit subroutine (either SPLINE, SPLINT, SPLISL, or SPINSL) will print this mes- 
sage, and the program will terminate. Since spline curves are used so extensively in 
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the program, it is difficult to state the possible cause. However, the printout of the 
spline points should assist in pinpointing the cause of the error. 


SPLINE ERROR -- NUMBER OF SPLINE POINTS GIVEN IS LESS THAN TNO 


NUMBER OF POINTS = -3 

X ARRAV Y ARRAY 

0.00003 10.000 


At least two points must be given to determine a spline curve. If not, the spline fit sub- 
routine (SPLINE, SPLINT, SPLISL, or SPINSL) will print this message, and the pro- 
gram will terminate. The one printed spline point may give a clue as to the error. 


NUMERICAL EXAMPLE 


An example of an axial- flow compressor rotor is used to illustrate the use of the 
program and to show the type of results that can be obtained. Flow is analyzed on the 
midchannel stream surface of an axial-flow rotor designed with the computer program 
of reference 11. The design pressure ratio is 1.275, the inlet hub-tip radius ratio, 

0.5; the aspect ratio, 1.5; the tip solidity, 1.0; and the tip relative Mach number at 
the inlet, 0. 9. Although tip relative Mach number is near sonic, there are no locally 
supersonic regions on the midchannel stream surface. 

The printed output from this example is given in table I. The input for this exam- 
ple is item 1 in that table. The number of mesh points used was 861: 41 in the axial 
direction, and 21 in the radial direction. 

This rotor was run twice, with reduction factors (REDFAC) of 0. 99999 and 0.70. 
This was done in order to test the effect of reduction factor and also the accuracy of the 
velocity- gradient method in comparison with the finite- difference method. Since both 
cases have REDFAC <1.0, both the subsonic (finite difference) and the transonic 
(velocity gradient) solutions will be obtained in each of them. (If REDFAC = 1.0, only 
the finite- difference solution is obtained.) 

The first case (REDFAC = 0. 99999) permits a comparison to be made between the 
approximate velocity- gradient method (used to obtain transonic solutions) and the more 
exact finite- difference, stream-function solution (used for subsonic solutions). Since 
the reduction factor is 0. 99999 (< 1 . 0) , both solutions will be obtained, but the answers 
should be very close if the methods compare well. 
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(c) Tip section. 

Figure 15. - Midchannel velocities for axial-flow compressor example. 


The second case (REDFAC = 0. 70) permits a comparison of reduction- factor effect. 
Since the finite- difference solution is now obtained at only 70 percent of the mass flow 
and rotational speed (appendix E), the results for the full- mass-flow, velocity-gradient 
solution will be less accurate than when a higher reduction factor is used. 

Figure 15 shows midchannel flow surface velocities for both reduction factors and 
both types of solution (the finite difference and the velocity gradient). Results are pre- 
sented for three values of stream function (0. 0, 0. 5, and 1.0) at hub, mean, and tip. 

On all parts of this figure, the solid line represents velocities from the final finite- 
difference iteration of the REDFAC = 0. 99999 solution. This solution is mathematically 
most accurate. The results presented here show a much better comparison of velocity- 
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MERIDL (REDFAC - 0.99999) 



la) Hub section. (b) Tip section. 

Figure 16. - Blade surface velocities for axial-flow compressor example. 


gradient to finite- difference results than that presented previously for the MERIDL pro- 
gram in reference 4. The two different velocity- gradient solutions are likewise now in 
much better agreement. 

Figure 16 shows a comparison of blade surface velocities estimated from the 
MERIDL program (appendix G) and calculated from the TSONIC program (ref. 3) . The 
TSONIC velocities are more accurate for the design blade shape. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 28, 1976, 

505-04. 
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APPENDIX A 


GOVERNING EQUATIONS 


The cylindrical coordinate system is shown in figure 4. However, the solution is 
obtained on an orthogonal mesh, as illustrated in figure 3. The variables s and t are 
used to denote distance along the streamwise and normal orthogpnals, respectively. 

For subsonic solution, the stream function is used. In appendix B we obtain the 
stream-function differential equation derived from the assumptions given in the section 
METHOD OF ANALYSIS: 


3 2 u 9^u _ 9u /sin <p + 1_ 9B + 1. 9£ _ 9<A _ 9u / cos (p + 1 _ 9B + J_ + &<£. 
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Equations (Al) to (A4) are derived from the hub-shroud momentum equation in 
appendix B. Note that all the partials in equations (Al) to (A 4) are on the stream sur- 
face except for 9p/90 in equation (A 4) , which is at constant z and r. 

The derivatives of the stream function satisfy the equations 

„ rBpW. 

9u = _ K t ( i 

9s w 
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(A 6) 


8u _ rB P W s 
at w 

For the final transonic solution the velocity-gradient method is used. The stream- 
line curvatures and flow angles needed for the velocity-gradient equation are obtained 
first from a reduced-flow, subsonic-stream-function solution. 

The velocity- gradient equation is 


dW = 


a w +b +— + d cos /3 ] dt +— + Wf 


W 


W 


(A7) 


The coefficients, a, b, c, and d, are given by different expressions in the blade re- 
gion, in the upstream region, and in the downstream region. These coefficients are 
given as follows: 


Blade- region coefficients: 


= £211 cos < a - <Pl - ^Li cos .£ + sin g sin 0 cos 0 £2 
r r 3t 


-\ 


dW 


dW, 


b = cos (3 — sin (a - cp) - 2 oj sin/3 cos <p + r cos /3( 2 + 2co sin al — V (A8) 


dm 


dm 


at 


c — 0 


d = 0 


Upstream- region coefficients: 
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Downstream- region coefficients: 



Finally, in all three regions, we have 


DT f? 

e = C T! - oj dA - C dT” +iii— dp” 

pi p n f * 


dT” R dp” 
2T ” 2CpP M 


(All) 
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Equations (A7) to (All) are derived in appendix C. Equation (A7) is solved as an 
initial-value problem, where the initial value of W is specified at the hub for any given 
t-line running from hub to tip. By finding several solutions for varying values of W at 
the hub, a solution satisfying continuity will be found; that is, the solution will satisfy 



pWrB cos^ - q>) cos 0 dt = w 


(A 12) 


When equation (A7) has been solved, subject to satisfying equation (A12), for every hub- 
to-tip mesh line in the region, the entire velocity distribution is obtained. 

The solution obtained by either the finite-difference or velocity- gradient method is 
for the midchannel surface between the blades. Of greater interest are the blade sur- 
face velocities. These can be estimated since the blade loading depends on the rate of 
change of whirl. By assuming a linear variation of velocity between blade surfaces, 
velocities can be calculated. In appendix G the following equation for calculating blade 
surface velocities is derived: 


B d(rV 0> 

W 7 = W . , - - COS 0 - 

l mid 2 P dm 


R d < rV fl> 

W tr = W mid + - cos 0 

tr miQ 2 dm 


(A13) 
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APPENDIX B 


DERIVATION OF STREAM- FUNCTION EQUATION 


Wu derives the following radial- and axial- momentum equations (eq. (96), ref. 1) 
on a meridional stream surface: 


r 9r 
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(Bl) 


The partial derivatives here all refer to the rate of change on the midchannel stream 
surface. Wu used a bold derivative symbol to indicate this, but it is not necessary here 
since we consider the equation only on the stream surface. The first equation of (Bl) is 
the radial momentum equation and the second is the axial momentum equation. 

Our flow solution is obtained by using finite- difference methods on an orthogonal 
mesh. The orthogonal mesh distances are s in the throughflow direction and t in the 
hub- tip direction. Figure 3 illustrates the s and t distances and also shows the angle 
cp between the s-direction and the axial direction. We want the momentum equation in 
the t-direction. This t-momentum equation will be used to derive the desired stream- 
function equation. The t-momentum equation is obtained by adding the t- components of 
the radial and the axial momentum; that is, by multiplying the first equation of (Bl) by 
cos cp and subtracting the second equation of (Bl) multiplied by sin cp. Also, the chain 
rule for a partial derivative is 


a a . a 

— = cos sin cp — 

at 8r 9z 


Therefore, the momentum equation in the t-direction can be written 
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(B2) 
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where is the t-component of the vector F, which is given by equation (97) of ref- 
erence 1 as 


F 


t 


n, r, 

t_ 9p 

n^rp 90 


(B3) 


Here n^ and n^ are the t- and 0-components of the unit normal to the stream surface. 
Note that 9p/90 is in the blade- to- blade direction at constant r and z and not on the 
stream surface as theother partials are. In the t-direction, the stream surface has an 
angle to the meridional plane. The tangent of this angle is r(90/9t). Hence, the com- 
ponents of the normal to the stream surface must satisfy 


n 


0 _ 


n t 




\9t, 


(B4) 


We can substitute this into equation (B3) to get 


F =— I ^2 
1 9t p 90 


(B5) 


It is desired to express the t- momentum equation (B2) entirely in s- and 
t-components and their derivatives. First, we will use 


W = W cos cp + W sin <p 
s z r 


(B6) 


This leaves only 


9W 9W 

r z 

9z 9r 


in equation (B2) to be expressed in s- and t- components. For this we use 

sin m 
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W = W, cos cp + W sin cp 
r t G ^ 
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W = W cos cp - W, sin cp 

Z S L 


(B7) 
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.nd the chain rule (noting that ds/dr = sin <p, etc.) 


r\ 3 0 

— = sin <p — + cos (p — 

at 
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(B8) 




3y using equations (B7) and (B8) we can express, after combining and simplifying, 
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The desired ^momentum equation in s and t could now be given by substituting equa- 
tions (B6) and (B9) into equation (B2) . 

It is desired to obtain the stream-function equation expressing the velocity compo- 
nents as partials of a stream function. This will reduce the number of unknowns in the 
momentum equation from two (W g and W t ) to one (the stream function) But it changes 
the equation to a second-order partial differential equation. A normalized stream func- 
tion u will be used. That is, u is zero at the hub and 1 at the shroud. The value of 
the stream function u can be obtained by integrating along a vertical mesh line. 




rB dt 
s 


By differentiating this expression, 

w 

rw S 


Similarly, the stream function u can be obtained by integrating along a vertical mesh 
line to any given mesh point s ][ , and then integrating along a horizontal mesh line: 



We can differentiate this expression with respect to s to obtain 
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£h = _i2£ w 

9s w 1 


Hence, the derivatives of the normalized stream function in the s- and t-directions are 
related to the velocity components by 


W 

9s w 

= w 

9t w b 


"A 


> 


J 


(BIO) 


The quantity B is the stream- surface thickness corresponding to the overall mass 
flow w. However, the variation in local stream sheet thickness is not known, so that 
the overall blade- to-blade spacing is used, and w is the mass flow through one blade 
channel. B is a function of z and r, so that we can use 


B(z, r) = 0 tr (z, r) - 0^(z,r) 


(Bll) 


where 0£ r and 0 ^ are the 0-coordinates of the trailing and leading blade surfaces, re- 
spectively . 

By solving equation (BIO) for and W g and differentiating, we obtain 
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Equations (BIO) and (B12) are now substituted into equation (B9) to eliminate W and 
W t in favor of the stream function u: 
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Now substitute equations (B6) and (B13) into equation (B2) and rearrange to obtain 
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The entropy change dS can be calculated from 

dS=C dT.Rife 
P T p 

Since the entropy of a particle at actual flow velocity is the same as the entropy of that 
particle at stagnation conditions (either absolute or relative), we can also use 

dS = C -£2111 - R d P - (B15) 

P p M 

This takes less calculation than using static values. 

The static temperature T is calculated from T = T" - W 2 /2C p . Using this with 
equation (B1 5), since I=C p T!-cu\, yields 

a^ + T as = a(C p T! - co\ - c p T”) RT „ 9p ,. + w^/r_ a^ _ _i_ 9 tA 

"at 9t at p” at 2 lc pP " at t" at J 


81 



But 
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Now, substitute equation (B16) into equation (B14) to obtain 
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where 
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and F^. is given by equation (B5) . 

In the program, sin <p and cos <p are stored, so that it is more convenient to cal- 
culate the partials of <p in equation (B17) by 


82 



a c£ _ 
as 

i 

a (sin cp) _ 

i 

a (cos cp) 

cos 

cp as 

sin co 

as 

II 

<t|£ 

1 

a (sin cp) _ 

1 

a (cos cp) 

cos 

cp 9t 

sin cp 

at 




(B20) 
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Either the middle term or the right term is used, depending on whether cos cp or 
sin cp is near zero. 

In equation (B 17), F t and 8(rV 0 )/8t must be calculated differently outside the 

blade and within the blade. Within the blade, 8(rV fl )/8t is calculated by a finite- 
difference approximation using the previous iteration. For the first iteration, the par- 
tial are assumed to be zero. Outside the blade, the partials are calculated from the 
input whirl distribution. From the input, the whirl is specified (or can be estimated) 
as a function of the stream function u. Hence, upstream of the blade, 
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Similarly, downstream of the blade, 
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The F t is caused by the hub-shroud pressure gradient induced by the blade lean. 
We can calculate F fc from equation (B5) . The blade-to-blade pressure gradient dp/dO 
can be calculated from the blade loading by assuming constant entropy blade to blade, 

so that 


dp = pC p dT 

Using the fact that T = T u - W^/2C p , we get 

9E = _ p w 

a Q dd 


since T” is constant from blade to blade. 

Now, substitute this in equation (B5) to get 


(B23) 


(B24) 
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The blade- to- blade velocity gradient is calculated from d(rV^)/dm, as explained in ap- 
pendix G, by using equation (G2) . Outside the blade, there is no blade- to- blade pres- 
sure gradient. Hence, 


F t = ° 


(B2G) 


outside the blade. 



APPENDIX C 


DERIVATION OF VELOCITY-GRADIENT EQUATION 


The general velocity-gradient equation is an expression for the value of the direc- 
tional derivative of the relative velocity. The velocity-gradient equation is derived 
from Newton's three-dimensional force equation. In this program the velocity- gradient 
equation is used to determine the approximate velocity distribution when there is locally 
supersonic velocity on the meridional midchannel stream surface. 

The general velocity-gradient equation can be written 
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Equation (Cl) is the same as equations (B13) and (B14) of reference 2 but with the addi- 
tional term T(dS/dq), to allow for variation in entropy. The entropy variation is known 
from the input specification of relative total-pressure-loss distribution. Entropy varia- 
tion along streamlines for adiabatic flow is possible only with viscous forces. However, 
these forces are primarily in the streamwise direction so that it is reasonable and con- 
sistent to neglect viscous forces along the hub-shroud mesh lines. 

Outside the blade row, c = 0. This follows from the fact that d(rV^)/ dm = 0 in this 
region. Using this identity and the relation = V^ - cor, we can derive the expres- 
sion 
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sin a 


dW 6> W 0 + 2wr 
dm r 

When this is substituted into the expression for c, we find that c = 0. 

The general velocity-gradient equations (Cl) and (C2) will be applied along vertical 
mesh lines in the t-direction, so that q = t. (See appendix B and fig. 3 for a descrip- 
tion of the orthogonal mesh.) Then 
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Note that, by using equation (C3), 
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Using equations (C3) and (C4) in equations (Cl) and (C2), we can get 
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= 1- i i-i- + cos (3 — sm(tt - cp) - 

dt r dm 

c 


W sin B , 0 . „ I 

1. + 2cu sin /3 


. d0 , 1 dT i dX „ dS ' 

dt w \ P dt dt dt 


(C5) 


We can express T dS in terms of the relative velocity W and relative stagnation 
conditions. We use 
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T - T „ vr 

2C. 


(C6) 


By using equations (B15) and (C6) , we get 


T dS = C_ dT” - dp” +— (— d P" - C 


dT” 


C p \ 2p” P 2T" 


(C7) 


Now we can write the velocity- gradient equation (C5) as 


dW = (aW + b)dt + — + Wf 
W 


(C8) 


where 


.2 


= cos p cosjq- 2) _ sin p cos £ + sin sin/3 cos p 


dO 

dt 


dW 


'dW 


b = cos /3 — sin(o - <p) - 2oj sin /3 cos 0 + r cos (3 | 

dm 


e 




,0 - 1 dO 

+ 2co sm o I — 

dm / dt 


>(C 9) 


e = C dT! - to dX- C dT” + BUZ dp” 
Pi P P M 


f - dT” _ R dp” 

2T” 2C p” 
P 


Equations (C8) and (C9) can be used directly within the blade row. However, out- 
side the blade row. the angle (3 is not known directly so that we desire to eliminate (3 
in favor of X, r, and W. This can easily be done since 


W 0 = -- wr 
r 


and 
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Then cos (3 is obtained from 


sin (3 



cos ^/3 = 1 - sirPp 

Also, c = 0 in equation (C2) outside the blade as was shown previously. With these 
substitutions, equation (C5) becomes 


dW 


( aW + — + d cos 
\ W 


P dt 


e 

W 


+ Wf 


(CIO) 


where 


a 


_ co s (q - t£) 


~\ 


c 



/ v . A. + cor 

cos(o / -</?) + 

r 


cos ip 


dW 

d = — sin(Q' - ip) 

dm 


e = C dT! - u> dA. - C dT" dp" 

Pi P 


> (Ctl) 


f _dT" R dp 1 * 

2T" 2C p" 
P 


cos (3 = 
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Downstream of the blade, outlet whirl should be used. Hence, in equation (Cll), 
/ r V^ is used instead of Equations (A7) to (All) are obtained directly from equa- 
tions (C8) to (Cll). 
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APPENDIX D 


LOSS CORRECTIONS 


An approximate loss correction is made by reducing the relative ideal stagnation 
pressure. This loss in stagnation pressure is specified by the input at the downstream 
input station by specifying either the loss fraction or the stagnation pressure. These 
quantities may vary from hub to tip. If the outlet stagnation pressure is given as input, 
the fractional loss of stagnation pressure is calculated. So, in all cases, the fractional 
loss of stagnation pressure at the outlet input station is known from hub to tip. 

For the case where the downstream stagnation pressure is given as input, the frac- 
tional loss of stagnation pressure is calculated. First, the outlet total temperature T' 
is calculated from the change in whirl along a streamline, 

“ R ,v <') 0 - A 1 

T' = T! + — t - — (Dl) 

o 1 c 


Then the ratio of actual to ideal stagnation pressure is calculated from 


P T 

o 



ideal 



(D2) 


Within the blade, the loss is distributed linearly or as specified bv the input, from zero 
at the leading edge to the fraction given or calculated at the blade trailing edge. The 
loss fraction is the same whether it is expressed in terms of relative or absolute stag- 
nation pressures. That is, 


,P T . T ., ^ .P L 

^ideal ^ideal 


(D3) 


When this ratio is known, the density can be calculated as a function of the velocity W. 
as follows: 
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.1/ (y-i) 


(D4) 



where 

T _ T! _ W 2 +2co\ - (cor) 2 (D5) 

1 2 C 

P 

Equation (D4) applies the loss correction by calculating the density from the ideal den- 
sity. Equation (D4) is derived as equation (B6) in reference 13. Equation (D5) is equa- 
tion (3) of reference 2. The value of p"/p£j eal varies at each point of the region and 
is recalculated with each iteration. The value of p is used in checking continuity by 
means of equation (A 12). 
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APPENDIX E 


DEFINING REDUCED- MASS- FLOW PROBLEM 


When the midchannel stream surface has locally supersonic flow, the solution can- 
not be obtained directly by solving the stream-function equation (eq. (Al)). However, 
an approximate solution can be obtained by getting a reduced-flow solution and extending 
this to the full- flow solution by the velocity-gradient method. This full-flow solution 
depends strongly on the reduced-flow solution, so it is important to establish the condi- 
tions that will give the most suitable stream- function solution. With a reduced mass 
flow, if both the inlet and outlet whirl and the rotational speed are reduced in the same 
ratio as the mass flow, a similar flow will result. That is, the reduced-flow stream- 
lines will have approximately the same angles and curvatures as the full- flow stream- 
lines throughout the flow field. 

Another consideration is the boundary conditions at the upstream and downstream 
boundaries. Again, the flow distribution will be similar to the reduced mass flow if the 
mass flow, inlet and outlet whirl, and rotational speed are all reduced by the same 
ratio . 

In summary, then, for the reduced- mass-flow solution, the mass flow, inlet and 
outlet whirl, and rotational speed are all reduced in the same ratio. This reduction 
factor is specified in the input by REDFAC, as explained in the section INPUT. 
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APPENDIX F 


INCIDENCE AND DEVIATION CORRECTIONS 

The solution region is divided into three subregions: upstream, within the blade, 
and downstream. These three regions should match each other; that is, there should 
be no difference in the tangential momentum of a particle on either side of the lines at 
the leading and trailing edges. If there were zero incidence and deviation angles, this 
condition would be satisfied automatically by using the actual blade shape in the blade 
region. Since there are almost always nonzero incidence and deviation angles, the 
mean flow will not follow the blade surface near the leading and trailing edges. Be- 
cause of this, an empirical correction is made to the blade shape near the leading and 

trailing edges. 

The blade- shape correction is made for a distance from the leading or trailing 
edge. This distance varies between l/6 and l/2 of the blade chord, depending on solid- 
ity. The solidity is defined to be the ratio of the true blade chord along a streamline to 
the blade- to- blade distance. When there is a change in radius, the solidity differs at 
the leading and trailing edges. For calculation purposes in the program the true blade 
chord along an s-coordinate line is used instead of that along a streamline. For a solid- 
ity of 2 or more, the distance is l/6 of the chord; and for a solidity of l/2 or less, the 
distance is l/2 of the blade chord 0 The distance varies linearly with solidity between 
these limits. The correction is made to the flow angle so that tangential momentum is 
continuous at the blade leading or trailing edge, and then the angle correction is lin- 
early decreased over the prescribed distance within the blade. This correction must be 
made to the flow angle j3. However, the blade shape in the stream-f unction equation is 
expressed by the partial derivatives, 30/ds and 30/dt. Therefore, with the corrected 
, the appropriate correction is made to 30/ds . It is assumed that the s-coordinate is 
close to a streamline. Therefore, only 30/ 3s is changed, and 30/3t is not changed. 

The free-stream flow angle is not known in advance at the blade leading edge. 
Therefore, an iterative procedure is used to make the blade- shape correction. On the 
first iteration, no blade-shape correction is made. After the first iteration, an approx- 
imate solution is obtained, and the flow angle /3 is calculated throughout the region. 
From the requirement of continuous tangential momentum, 

( w *) fs = W bf 


From continuity, 


C-l' 
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where NBL is the number of blades. Hence, since tan/3 = W^/ W , 


tan /3 bf = tan /3 fg 


B x NBL p bf 
27r p fs 


(FI) 


Equation (FI) is used to calculate f}^. The difference /3 b j - /3 b (see fig. 12) is then 
varied linearly within the blade for the specified distance. We can express /3 in terms 
of the partials of 0 by 


tan/3 bf 




cos(ci! - (p) +— sin(o; - 

at 


<p) 


When this is solved for (80/as) b £, we obtain 

tan ' J bf 39 . , 

/M _ r 3t (FS 

\as/ bf cos {a - (p) 

Equation (F2) is used to calculate (aO/as) b £ from the linearly varying /3 at each mesh 
point. The calculations for these blade-shape corrections are done by subroutine 
INDEV for each iteration. When REDFAC is less than 1, the final transonic calcula- 
tion uses the same method to correct /3 within the blade (subroutine LINDV), but 
a 0/as is not needed. 
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APPENDIX G 


BLADE SURFACE VELOCITIES AND BLADE -TO- BLADE AVERAGE DENSITIES 


The blade surface velocities can be calculated once a meridional midsurface solu- 
tion is obtaihed. The blade- to- blade velocity gradient 3W/30 depends on dW^/dm. 

To obtain the desired relation, we use the velocity- gradient equations (Cl) and (C2). In 
this case, q = 9 and dr/d0 = dz/d0 = 0. Also, it is assumed that T!, and the en- 
tropy S are constant in the blade- to- blade direction. With this, equations (Cl) and (C2) 
become 


3W 

= W sin n sin jS cos j3 + r cos /3 

do 



+ sin O' 


(Gl) 


Using the fact that V 0 = W 0 + cor and W 0 = W sin /3 , we can rewrite equation (Gl) as 


3W 
3 9 


cos f3 


d(rV 0 ) 

dm 


(G2) 


Since only a midchannel solution is obtained, we assume that cos /3 and d(rVg)/ dm 
are constant from blade to blade. This means that W varies linearly from blade to 
blade. (This may be in considerable error near the leading or trailing edge. ) With this 
assumption, then 


W, + w 

W . , =— i — 

mid o 


(G3) 


Integrating (G2) from blade to blade (from 9 = 6L to 9 = 0^ and using (G3), we obtain 


W. = W . , + 
tr mid 


B cos i 3 d(rV 3 ) 




dm 


W , = W . , 
l mid 


B cos B d(rV ^ 


dm 


(G4) 


Note that equation (G4) is very close to that developed by Stanitz (eqs. (16) and (17), 
ref. 16) . 
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Equation (G4) is used in subroutine BLDVEL to obtain blade surface velocities from 
the midchannel meridional solution. 

It is desirable to consider the blade- to-blade variation in density to satisfy continu- 
ity for the blade passage. Equation (G4) gives information we need to do this. Since the 
midchannel solution is considered to be representative of all meridional-plane stream 
surfaces, we can consider the solution to be based on average blade— to— blade conditions. 
This means that in equation (BIO) average blade- to-blade velocity and density should be 
used. Equation (BIO) becomes 


— = £2 

3s w 


K) 

\ vav 




9u _ _rB 
9t w 


K) 

V /av 


(G5) 


If we use Simpson r s male to calculate the average value, we obtain 


/ W ^ + 4p mid W s, mid + p tr^s, tr 

' s 4v 6 


(G6) 


The velocity component W g can be expressed in terms of the velocity and flow angles 

by 


W = W cos ft cos (o - tp) 


(G7) 


This relation holds also on the blade surfaces, since the flow angles are assumed to be 
constant in the blade- to-blade (0) direction, so that 


By using equations (G3) , 


W g ^ = \V L cos (3 cos (o - cp) 


(G8) 


W 


tr = W t cos (3 cos (o - (p) 


(G4) , and (G8) in equation (G6), we can express (pW ) 

V s 4v 


as 


(pw s ) av = p a w s 


Pi Ptr 

+ — cos 13 cos (O' - <p)(Wj - W tr ) 


(G9) 
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where the subscript mid is omitted from W , and 

S 

+ 4 P m id + P tr 


WTien a solution of equation (Al) is obtained, all quantities in equation (G9) can be esti- 
mated from the previous iteration, except for W . So we solve for W in equa- 
tion (G9) and use equation (G5) to obtain 


W 

s 


9u 
w — 

at 

rBp 
' av 


(p£ ~ P tr ) cos (3 cos (O’ - (ft) (W^ - W ) 

12p 

p av 


(CIO) 


ITiis is the equation used in subroutine NEWRHO to calculate W from 9u/0t. On the 
first iteration, p ay = p '^ ea j is used, and the second term is omitted. After this, the 

values of p^^, p£ r , p^ , W^, W^, cos/3, and cos(» - cp) from the previous iteration 
are used. 

In a parallel manner, we can derive the equation for W^: 


W t = - 


rBp 


Ou 

9s <P Z " p tr ) COS|S sin(Q " ^ (W /. “ W t,-> 

12p, 


av 


av 


(Gil) 
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APPENDIX H 


SYMBOLS 

a coefficient in velocity- gradient equations (A7) and (Cl) 

B tangential space between blades , rad 

b coefficient in velocity- gradient equations (A7) and (Cl) 

Cp specific heat at constant pressure, J/ (kg)(K) 
c coefficient in velocity- gradient equations (A7) and (Cl) 

d coefficient in velocity- gradient equation (A7) 

e coefficient in velocity- gradient equation (A7) 

F vector normal to midchannel stream surface and proportional to tangential pres- 
sure gradient, N/kg 

f coefficient in velocity- gradient equation (A7) 

2 2 

I rothalpy, C^T! - coX, meters /sec 

i incidence angle, rad 

m meridional streamline distance, meters 

n unit vector normal to midchannel stream surface 

p pressure, N/meter^ 

q distance along an arbitrary space curve, meters 

R gas constant, j/(kg)(K) 

r radius from axis of rotation, meters 

r radius of curvature of meridional streamline, meters 

c 

S entropy, J/(kg)(K) 

s distance along orthogonal mesh lines in throughflow direction, meters 

T temperature. K 

t distance along orthogonal mesh lines in direction across flow 

u normalized stream function 

V absolute fluid velocity, meters/sec 

W fluid velocity relative to blade, meters/sec 
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w mass flow, kg/ sec 

z axial coordinate, meters 

a angle between meridional streamline and axis of rotation, rad; see fig. 4 

/3 angle between relative velocity vector and meridional plane, rad; see fig. 4 

y specific-heat ratio 

£ coefficient in stream- function equation, defined in eq. (A3) 

0 relative angular coordinate, rad; see fig. 4 

A prerotation, ^rV^j , meters^/sec 

£ coefficient in stream- function equation, defined in eq. (A2) 

p density, kg/meter 

cp angle between s-distance line and axis of rotation, rad; see fig. 3 
ip stream function, kg/ sec 

rotational speed, rad/ sec; see fig. 4 
Subscripts: 

av average blade- to- blade 

b blade 

bf blade flow 

cr critical 

fs free stream 

h hub 

1 inlet 

l blade surface facing direction of positive rotation 

m component in direction of meridional streamline 

mid midchannel blade- to- blade 

o outlet 

r component in radial direction 

s component in s- direction 

t tip, or component in t- direction 

tr blade surface facing direction of negative rotation 
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z component in axial direction 

6 component in tangential direction 

Superscripts: 

* absolute stagnation condition 

" relative stagnation condition 
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